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COMPUTER-AIDED  DESIGN  OF  INTEGRATED  CIRCUIT 


FABRICATION  PROCESSES  FOR  VLSI  DEVICES 

SUMMARY 

The  work  described  here  builds  upon  an  established  program  which  over 
the  past  five  years  has  resulted  in  the  development  of  three  generations  of 
computer-aided  process  modeling  programs  -  SUPREM  I,  II  and  III.  The  second 
generation  version  -  SUPREM  II  -  is  currently  in  use  at  more  than  150 
industrial  firms  and  university  research  groups  and  is  generally  regarded  as 
the  preminent  tool  of  its  type  available  today.  A  significant  number  of 
published  papers  from  groups  all  over  the  world  reference  this  computer 
program,  und  in  fact,  many  papers  include  simulated  results  obtained  from 
SUPREM.  It  has  found  widespread  application  in  the  design  of  state-of-the- 
art  MOS  and  bipolar  fabrication  processes. 

The  principal  highlight  of  the  program  during  the  present  reporting 
period  has  been  the  release  of  SUPREM  III.  This  is  a  completely  new,  greatly 
improved  version  of  the  program  with  substantially  improved  capabilities. 
SUPREM  III  can  simulate  structures  with  up  to  five  material  layers.  Built  in 
models  for  Si,  Si02,  Polysilicon  and  Si 3N are  included  and  other  materials 
can  easily  be  added  with  user  defined  parameters.  A  large  number  of  new 
and/or  improved  physical  models  for  oxidation,  diffusion,  epitaxy  and  ion 


implantation  are  included  in  SUPREM  III,  these  have,  for  the  most  part,  come 
from  experimental  and  theoretical  work  under  this  program  over  the  past  five 
years.  These  new  models  should  greatly  expand  the  range  of  applicability  and 
the  accuracy  of  the  program. 

At  the  present  time  approximately  25  copies  of  the  new  program  have  been 
licensed  from  the  University  by  U.  S.  semiconductor  manufacturers  and  defense 
contractors.  We  expect  that  over  the  lifetime  of  the  program  the  distribution 
will  eventually  exceed  that  of  SUPREM  II.  Initial  feedback  from  industry 
users  is  beginning  to  come  in  to  us,  and  it  is  clear  already  that  the  new 
program  capabilities  will  make  SUPREM  III  a  much  more  valuable  tool  than  its 
predecessors.  We  expect  over  the  next  two  years  to  release  updated  versions 
of  the  program  as  new  materials  models  become  available.  Much  of  the  work 
which  is  aimed  at  this  objective  is  described  in  the  remainder  of  this  report. 

The  specific  goals  of  this  research  have  been  and  continue  to  be  the 
formulation  of  basic  physical  and  mathematical  models  of  integrated  circuit 
fabrication  processes  which  can  accurately  predict  the  structure  of  small 
geometry  VLSI  devices  (and  other  arbitrary  device  structures)  which  will 
result  from  a  given  fabrication  sequence,  and  to  implement  these  models  in  a 
comprehensive  process  modeling  computer  program.  Such  a  program,  it  was 
believed  at  the  outset  of  this  work,  would  significantly  reduce  costly  and 
time  consuming  iterative,  experimental  approaches  to  developing  or  optimizing 
silicon  technologies.  This  has,  in  fact,  occurred.  SUPREM  III  represents  a 
major  step  along  the  path  toward  improved  simulation  tools. 

It  is  clear  that  the  historical  trends  in  integrated  circuit  technology 
over  the  past  two  decades  towards  increased  complexity  and  smaller  active 
device  dimensions  will  continue.  Modern  device  structures  employ  lateral 
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device  dimensions  on  the  order  of  2-3  p  and  vertical  dimensions  well  below 
1  pm.  There  are  no  basic  physical  mechanisms  which  will  prevent  a  reduction 
in  each  of  these  dimensions  by  an  additional  order  of  magnitude  over  the  next 
two  decades.  Practically  accomplishing  this  continued  scaling,  however, 
depends  upon  our  ability  to  physically  understand  and  quantitatively  model  the 
fabrication  techniques  which  will  be  used  in  the  construction  of  such  devices. 

For  devices  with  relatively  large  geometries  (>5  p)  and  loose  processing 
tolerances,  relatively  simple  models  suffice  for  prediction  of  vertical  device 
structures  resulting  from  a  given  fabrication  sequence.  As  device  dimensions 
shrink,  however,  it  becomes  essential  to  employ  more  robust  process  models 
and  to  consider  the  interaction  both  laterally  and  vertically  of  various 
processing  steps,  if  accurate  simulation  of  structures  is  to  be  obtained. 

This  is  important  even  with  today's  2-3  pm  device  geometries;  it  will  become 
essential  for  smaller  devices. 

Large  geometry  devices  can  be  successfully  modeled  as  one-dimensional 
structures.  This  is  true  for  both  process  models  and  electrical  models. 
Devices  with  lateral  dimensions  below  a  few  microns,  however,  require 
two-dimensional  models  for  accurate  simulation.  This  need  has  stimulated  a 
large  body  of  work  in  recent  years  on  two-dimensional  electrical  models  of 
device  current-voltage  characteristics.  One  example  of  this  type  of  work  is 
the  2D  program  GEMINI,  developed  under  this  program.  Work  of  this  type  has 
resulted  in  remarkable  advances  in  our  understanding  of  small  geometry  device 
physics. 

Progress  has  not  been  as  rapid,  however,  in  two-dimensional  process 
modeling.  While  some  basic  2D  process  modeling  programs  have  been  developed 
at  Stanford  (SUPRA  for  example)  and  elsewhere,  these  programs  do  not 


incorporate  robust  2D  kinetic  models  which  can  accurately  predict  doping 
profiles  and  device  geometries  under  a  wide  range  of  conditions.  This  is  a 
direct  result  of  our  need  for  improved  physical  models  of  oxidation,  ion 
implantation,  diffusion,  and  CVD.  It  is  quite  clear  that  these  processes  are 
not  one  dimensional.  Recent  experimental  evidence  has  clearly  indicated  that 
oxidation  or  impurity  diffusion  in  a  localized  region  of  a  silicon  substrate 
can  substantially  affect  oxidation  or  diffusion  rates  in  laterally  or 
vertically  adjacent  regions  of  the  substrate.  There  is  no  clear  agreement  at 
the  present  time  on  the  basic  physical  mechanisms  responsible  for  such 
results,  although  much  has  been  learned  in  the  past  several  years.  It  is 
clear,  however,  that  we  must  quantitatively  understand  such  phenomena  if  we 
are  to  accurately  model  small  device  structures. 

A  specific  goal  of  this  program  is  to  understand  and  model  these 
two-dimensional  effects.  We  have  made  substantial  progress  in  this  regard  in 
the  past  six  months.  It  appears  now  that  the  basic  physical  phenomena 
underlying  these  interactions  are  the  roles  of  points  defects--silicon 
vacancies  and  interstitials--in  impurity  diffusion,  thermal  oxidation,  and 
other  processes.  The  generation  and  consurrption  of  these  point  defects  during 
high  temperature  fabrication  steps  appear  to  be  the  unifying  physical  effects 
which  can  explain  many  of  the  phenomena  which  have  been  regarded  as  anomalous 
to  date.  We  have  used  such  models  to  quantitatively  understand  a  variety  of 
process  phenomena  and  have  incorporated  some  of  these  models  in  SUPREM  II. 

More  such  models  are  in  SUPREM  III.  The  unifying  role  of  these  mechanisms  is 
a  cornerstone  of  the  work  described  in  this  progress  report.  We  regard  such 
models  as  absolutely  essential  to  accurate  modeling  of  two-dimensional  effects 
in  small  devices. 
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Simply  stated,  the  long-range  objective  of  this  research  is  the 
development  of  a  2D  process  modeling  program  which  contains  physically  correct 
2D  kinetic  models.  Most  of  the  effort  under  this  program  is  aimed  at 
experimental  and  theoretical  work  to  uncover  the  basic  physical  mechanisms 
which  govern  2D  oxidation  kinetics,  2D  diffusion,  etc.  This  is  a  major 
problem  area  with  very  difficult  materials  problems.  It  will  become  clear  in 
the  technical  discussions  later  in  this  report  that  2D  process  modeling 
involves  much  more  than  solving  ID  equations  in  two  dimensions.  Basic 
physical  laws  and  understanding  are  lacking  at  present.  Effects  such  as 
lateral  OED,  oxidation  under  masking  Si3N4  layers,  lateral  and  vertical 
effects  of  high  dopant  concentrations  on  diffusion  coefficients,  and  a  host  of 
other  known  experimental  effects  cannot  be  explained  by  a  simple  extrapolation 
of  known  ID  physical  laws  to  2D  structures.  This  program  is  aimed  at 
generating  the  physical  understanding  needed  to  develop  2D  kinetic  models. 

An  essential  part  of  this  is  to  develop  models  for  bulk  point  defect 
(interstitial  and  vacancy)  generation,  recombination  and  diffusion,  since  it 
is  clear  that  the  local  concentrations  of  these  defects  determine  local 
diffusion  coefficients,  oxidation  rates,  etc.  In  fact,  an  alternative 
statement  of  the  overall  objective  of  this  program  would  be  to  develop 
techniques  for  calculating  local  (i.e.  time  and  position  dependent)  process 
parameters  suitable  for  process  simulation.  Such  process  parameters  will  of 
necessity  be  geometry  dependent  which  means  that  diffusion  coefficients  and 
oxidation  rate  constants,  for  example,  will  depend  on  the  presence  or  absence 
of  nearby  heavily  doped  regions  or  oxidizing  interfaces.  This  will  imply  a 
tight  coupling  between  surface  geometry  and  resulting  impurity  profiles  in 
small  devices. 
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In  this  context,  the  overall  objectives  of  this  program  are: 

(1)  To  develop  accurate  physical  models  for  the  basic  technologies  used 


in  the  fabrication  of  silicon  integrated  circuits— oxidation,  ion 
implanatation,  diffusion,  and  chemical  vapor  deposition  (CVD). 

(2)  To  develop  accurate  physical  models  for  the  interaction  of  these 
processes,  for  example,  the  enhancement  of  diffusion  coefficients  in 
an  oxidizing  ambient. 

(3)  To  develop  accurate  physical  models  based  upon  silicon  point  defects 
for  two-dimensional  process  phenomena  in  small  geometry  devi  .. 

(4)  To  develop  high  resolution  analytic  tools  for  the  experiment 
characterization  of  small  device  structures. 

(5)  To  implement  improved  process  models  in  updated  versions  of 

III  so  that  this  program  is  capable  of  accurately  simulating  small 
devices  and  multi-layer  structures. 

(6)  To  aim  toward  a  robust  2D  process  modeling  program  (SUPREM  IV)  which 
would  combine  outputs  from  this  program  with  numerical  techniques 
developed  under  other  support. 

(7)  To  disseminate  the  results  to  the  semiconductor  and  defense 
industries. 

The  remainder  of  this  report  describes  the  specific  research  activities 
during  the  six  month  period  covered  by  this  report.  Each  of  the  subsequent 
sections  includes  a  discussion  of  the  problems  under  investigation,  the 
results  to  date  and  plans  for  the  future.  The  unifying  themes  of  all  of  the 
work  to  be  described  are  the  development  of  new  or  improved  models  for  SUPREM 
III  and  the  investigation  of  2D  materials  phenomena  required  to  develop  SUPREM 
IV. 
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1.  OXIDATION  AND  POINT  DEFECTS 
S.  DUNHAM,  J.  PLUMMER 

The  existence  of  interactions  between  silicon  oxidation  and  the  diffusion 
of  impurities  (OED  and  ORD)  and  growth/shrinkage  of  stacking  faults  (OISF)  has 
been  well  established  [1-12].  Both  effects  are  attributed  to  the  creation  of 
interstitialcies  at  the  oxidizing  interface  in  order  to  relieve  the  stress 
created  there  [13-15].  Many  of  these  studies  have  examined  the  relationship 

dx 

between  oxidation  rate  o  and  interstitialcy  density  resulting  in  fits 

It- 

to  the  equation: 


Unfortunately,  the  value  of  n  must  be  calculated  empirically  and  reported 
values  have  ranged  from  less  than  0.3  [4]  to  more  than  1.0  [12].  Our  ongoing 
work  then,  is  an  attempt  to  physically  account  for  the  relationship  between 
oxidation  rate  and  point  defect  densities  and  thus  to  develop  quantitative 
models  for  OED,  ORD  and  other  phenomena  suitable  for  2D  process  modeling  in 
SUPREM  IV. 

It  has  been  observed  that  in  lightly-doped  silicon,  oxidation  is  more 
than  99%  complete  [16],  so  virtually  all  of  the  silicon  interstitialcies 
created  by  oxidation  must  diffuse  into  the  oxide  where  they  react  with 
incoming  oxygen  to  complete  the  process.  Therefore,  the  fluxes  within  the 
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oxide  (Figure  1)  determine  the  interstitialcy  density  at  the  interface.  We 
can  write  the  continuity  equations  for  oxygen  and  silicon  interstialcies  as: 


at 


'  kICICo 


(2) 


"  kICICo 


(3) 


In  steady  state,  as  for  relatively  long  oxidations,  the  time  derivatives  are 
zero  so  combining  (2)  and  (3)  and  integrating  twice  we  get: 

CQ  =  yCi  +  a  +  bx  (4) 

where  y  =  (Dj/OgO). 

Tracer  experiments  [17-19]  have  shown  that  oxidation  takes  place  very 
close  to  the  oxide-bulk  interface  so  there  is  no  flux  of  interstitials  through 
the  oxide.  Therefore,  we  can  conclude  that: 

(irl  x  =  o  =  0  (5) 

and  Cj(o)  =  0  (6) 

By  Henry's  Law:  CQ(0)  =  KPq  (7) 


so  combining  with  (4)  and  (6): 


a 


(8) 
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The  rate  of  oxide  growth  is  proportional  to  the  flux  of  oxygen  at  the  gas- 
oxide  interface: 


D  ldCo\  -  1  dXo 

o  1  liTj  x  -  o  q  ~3t~ 


(9) 


where  X0  is  oxide  thickness  and  q  is  the  molecular  volume  of  Si 0 2 - 
Combining  this  with  (4)  and  (5): 


b  =  -  (DQ) 


-1 


fr) 


(10) 


We  can  now  rewrite  equation  (4)  as: 


c0  -  rC,  ♦  kp02  -  (D^r1  /£), 


(ID 


For  conditions  where  the  partial  pressure  dependence  of  the  oxidation  rate  is 
linear  (wet  or  high  temperature  dry),  at  the  interface  the  rate  of  oxidation, 
which  is  equal  to  the  oxygen  flux,  is  proportional  to  the  density  of  oxygen: 


/dCo\ 

\dx  / 


x  =  x. 


■  kC  Xn 
000 


(12) 


Due  to  volume  constraints,  a  certain  percentage  (*)  of  interface  reactions 
must  result  in  point  defect  interactions.  Ho  and  Plummer's  results  [20,21] 
show  that  vacancy  mechanisms  are  negligible  in  intrinsin  silicon  so: 


0 


I 


x  =  x 


0 


(13) 
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Substituting  the  derivative  of  (11)  and  rearranging: 


(14) 


Combining  this  result  with  (12)  and  the  experimental  support  for  the 
Deal -Grove  model  [22]  that  says: 


X 


(15) 


then  solving  for  the  density  of  interstitialcies  at  the  oxide-bulk  interface 
we  get: 


(16) 

(17) 

(18) 


We  can  see  that  depending  on  the  magnitudes  of  and  K^,  which  will  change 

with  temperature  and  partial  pressure,  any  power  law  relationship  relating 
Cj  and  dx0/dt  between  0  and  1  can  be  approximated. 

Looking  first  at  stacking  fault  growth: 

t 

L  =  /  (KjCj  -  KyCy  -  KR)  dt  (19) 

0 

.  ★ 

Assuming  steady-state  recombination  of  point  defects,  then  Cy=  (Cy  C 
so  using  (16)  and  looking  only  at  parabolic  oxidations  where  4r-  * 
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we  can  integrate  to  get: 


Ignoring  the  final  term  which  should  be  small  since  Cj>>C^  ,  this  expression 

gives  a  stacking  fault  length  time  dependence  of  between  0.5  and  1.0  which 

1/2 

fits  well  with  reported  dependences  [3,7-11].  In  fact,  fits  to  l  *  A^t  +A2t 
are  virtually  indistinguishable  from  the  accepted  practice  of  fitting  to 
L  *  Atn  (Figure  2).  Since  the  constants  A1  and  A.,  will  change  with  oxidizing 
conditions,  we  can  see  how  different  apparent  power  relationships  should 
result. 

Most  impurities  diffuse  via  point  defect  mechanisms  which  in  the  most 
general  case  gives: 

Di  =  CID.I  +  CVDV 

As  for  stacking  faults,  except  for  impurities  that  diffuse  only  with  vacancies 
(ORD),  the  second  term  is  usually  negligible  so  calculating  the  time-averaged 
diffusion  constant  gives: 

dx 

<  °i  >  -  01K1  <  >  +  DjK2  (22) 

which  can  apear  as  any  power  dependence  between  0  and  1.  Most  experimental 
work  has  involved  changing  either  temperature  or  partial  pressure  so  no  direct 
conparison  with  experiment  has  been  done  to  this  point. 


In  conclusion,  this  work  gives  a  physical  basis  for  the  observed 
dependence  of  impurity  diffusion  and  stacking  fault  growth/regrowth  on 


oxidation  conditions.  There  appears  to  be  great  potential  in  unifying  the 
many  seemingly  contradictory  results  that  have  shown  up  in  this  area.  Such  a 
model  would  enable  prediction  of  the  effect  of  oxidation  on  the  resulting  bulk 
properties  and  perhaps  point  out  how  to  most  effectively  make  use  of  these 
interactions. 
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2.  THIN  OXIDE  GROWTH  KINETICS  IN  DRY  OXYGEN 
H.  MASSOUD,  C.  HO,  J.  PLUMMER 


Very  thin  layers  of  thermal  Si02  were  grown  in  dry  oxygen  in  the 
thickness  range  0  -  500  A  as  described  in  earlier  reports.  These  oxidations 
were  carried  out  on  single  crystal  silicon  of  different  orientations  (100), 
(111)  and  (110)  ,  for  different  phosphorous  substrate  doping  densities 

(1  x  1015cm“3),  (2  x  1020crrr3)  and  (3.2  x  1020cnr3)  ,  in  different 

oxygen  partial  pressures  in  argon  (1  atm,  0.1  atm  and  0.01  atm),  in  the 
temperature  range  800°C  to  1000°C.  The  growth  was  monitored  by  an  in-situ 

automated  high-temperature  ellipsometer  which  is  well-suited  for  data 

( 

collection  in  the  initial  stages  of  oxidation  where  the  growth  rate  is 
enhanced.  The  experimental  work  was  done  in  cooperation  with  IBM. 

The  data  were  analyzed  in  the  context  of  the  Deal -Grove  linear-parabolic 
model  where  the  oxidation  rate  is  expressed  as: 


dXox  B 

dt  '  2Xqx+A 


(1) 


where  B  and  (B/A)  are  the  parabolic  and  linear  oxidation  rate  constants 
respectively.  These  rate  constants  describe  the  oxidation  kinetics 
beyond  the  fast  initial  regime.  The  rate  constants  were  obtained  using 
the  optimum  X^  method.  The  oxidation  rate  enhancement  in  the  thin 
regime  can  be  expressed  as  an  excess  or  as  a  ratio  term: 


IT  "  Ix^  +  ^ox’^ 

cIT*  *  t1  +  W)) 


(2) 


(3) 
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where  Ag  and  Af  are  the  excess  or  ratio  enhancement  terms  which  vanish 
with  the  onset  of  the  linear  parabolic  kinetics.  For  all  the  orientations 
studied  in  100%  pure  oxygen,  it  was  found  that  these  excess  terms  could  be 
fitted  to  the  following  expressions: 


dXox 

"3t 


8 


Cl6 


ox/L, 


-  X 


C2e 


ox/L, 


(4) 


or 


dXox  . 
It — 


B 


★ 

Kle 


-  t/ 


"1 


* 

K2e 


-  t/ 


x2 ' 


(5) 


The  dependence  of  the  oxidation  excess  enhancement  term  a0  on  oxide-thickness 
is  defined  by  two  terms  exponentially  decaying  with  thickness.  Both  terms 
contribute  to  Ag  up  to  -  50A,  where  the  term  (Cj.Lj)  vanishes.  From  -  50 A  to 
to  the  onset  of  linear-parabolic  kinetics  (-  300A),  only  the  term  (C2,L2) 
is  responsible  for  the  rate,  enhancement:  L^is  of  the  order  of  7  -  17 A,  and  is 
a  weak  function  of  temperature.  L-,  is  of  the  order  of  60  -  80A,  is  tempera¬ 
ture  independent  and  is  only  a  function  of  substrate  orientation. 

The  oxidation  rate  enhancement  studied  as  a  ratio  term  was  best  fitted  to 
an  expression  like  that  of  equation  (5).  In  this  analysis  K*  and  K2  are 
orientation  dependent  constants  that  are  temperature  independent.  This 
analysis  is  most  consistent  with  the  assumption  that  the  concentration  of 
oxygen  species  in  the  oxide  during  the  initial  stages  of  the  oxidation  is 
enhanced  above  what  is  inferred  from  permeation  studies.  The  most  frequently 
quoted  values  for  the  solubility  of  oxygen  and  water  modules  in  Si02  are 
5  x  lO^cm"^  and  5  x  lO^cm"^  respectively.  This  corresponds  to  an 
average  spacing  between  the  oxidant  species  of  270A  and  27A  for  oxygen  and 
steam  respectively. 


-17- 


It  is  easily  seen  that  the  concept  of  an  average  bulk  solubility  is 
questionable  in  the  case  where  the  average  spacing  between  the  solute 
molecules  is  larger  than  the  thickness  of  the  film  considered.  This  is 
illustrated  in  Figure  1.  The  same  conclusion  is  arrived  at  when  we  consider 
the  validity  of  application  of  Henry's  Law  to  the  kinetics  of  oxidation  in  the 
initial  stages.  Henry's  Law  assumes  that  the  solubility  C  of  a  gas  in 
equilibrium  with  a  solid  is  related  to  the  gas  pressure  pQ  by  the  Henry's 
Law  constant  H: 

C  =  HP(J  (6) 

This  assumes  a  semi-infinite  solid.  This  assumption  is  needed  to  establish 
the  equilibrium  between  the  flux  of  gas  molecules  entering  the  solid  and  the 
flux  of  gas  molecules  leaving  the  solid.  In  equilibrium,  these  two  fluxes  are 
equal.  In  the  case  of  a  thin  layer  of  Si02  (as  shown  in  Figure  1),  it  is  seen 
that  in  the  case  of  dry  oxygen,  the  flux  of  oxygen  molecules  consumed  in  the 
oxidation  reaction  at  the  Si -SiO 2  inteface  may  directly  affect  the  equili¬ 
brium  established  at  the  outside  surface  of  the  Si 02  layer.  In  the  case  of 
steam,  the  effects  of  the  flux  of  water  molecules  consumed  at  the  Si -SiO 2 
interface  are  screened  by  the  higher  solubility  of  water  in  Si02.  This  is  in 
agreement  with  the  absence  of  a  fast  initial  oxidation  stage  in  the  case  of 
steam. 

Analysis  of  the  oxidation  data  at  lower  partial  pressures  of  oxygen  will 
help  elucidate  the  importance  of  the  role  oxygen  solubility  plays  in  the 
non-steady  state  period  present  in  the  case  of  oxidation  in  dry  oxygen.  It 
may  also  suggest  other  physical  mechanisms  that  play  a  role  in  fast  initial 
oxidation  as  well. 
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FIGURE  1:  Physical  arguments  which  may  explain  why  O2 
solubility  is  higher  than  the  equilibrium 
value  during  the  initial  oxidation  stage 
(see  text). 
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I.  ATMOSPHERIC  AND  HIGH  PRESSURE  OXIDATION  OF  POLYCRYSTALLINE  SILICON 

B.  CAIRNS,  L.  LIE 

High  pressure  oxidation  runs  for  polycrystall ine  silicon  kinetic  studies 
were  not  made  during  this  time  period  because  a  problem  of  accurate  tempera¬ 
ture  control  developed  in  the  high  pressure  oxidation  system  when  it  was 
converted  from  a  3- inch  to  4- inch  wafer  capability  for  other  internal 
Fairchild  programs.  A  new  tube  and  thermal  control  system  which  will  correct 
this  problem  have  been  ordered  and  high  pressure  runs  should  resume  in  the 
late  October  to  early  December  time  frame. 

3.1  Resistivity  of  Oxidized  Polycrystall ine  Silicon 

Measurements  of  polycrystalline  silicon  film  resistivities  following 
pyrogenic  steam  were  completed.  The  phosphorous  doped  polycrystalline  silicon 
fims  were  oxidized  in  1,  5  and  10  atm  pyrogenic  H20  at  750°,  800°,  and  850°C. 
Results  show  that  the  resistivity  is  dependent  on  oxidation  temperature  and  is 
fairly  constant  as  the  oxidation  proceeds  at  a  given  tenperature  (Figures  1 
and  2).  As  expected,  the  higher  oxidation  temperatures  yielded  lower  film 
resistivities.  In  all  cases,  we  did  not  oxidize  more  than  50%  of  the  poly-Si 
films.  Results  reported  by  Saraswat  of  Stanford  show  a  substantial  increase 
in  resistivity  when  more  than  60%  of  the  films  are  oxidized.  The  polycrystal¬ 
line  silicon  films  which  did  not  receive  a  pre-oxidation  anneal  show  a  very 
rapid  increase  in  the  film  resistivity  for  all  temperatures  and  pressures  to 
the  values  for  the  annealed  films.  This  is  not  surprising  for  oxidation 
temperatures  below  the  950°C  at  which  these  films  were  doped.  Reduced  dopant 
solubility  within  the  polycrystalline  silicon  grains  after  oxidation  could  be 
the  main  factor  for  this  rapid  increase  in  the  resistivity.  We  could  not 
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determine  a  time  constant  for  the  process,  but  it  is  complete  in  about  10 
minutes.  It  was  also  noted  that  pressure  effects  appear  negligible  on  the 
variation  of  film  resistivities  following  oxidations  up  to  10  atm  at  these 
temperatures.  The  slight  differences  between  1  atm  and  20  atm  seen  between 
Figures  1  and  2  are  within  the  range  of  data  scattering. 

3.2  Atmospheric  Oxidation  of  Polycrystalline  Silicon 

The  oxidations  of  undoped  polycrystalline  silicon  films  in  1  atm 
pyrogenic  steam  at  750°,  800°,  and  850°C  have  been  completed.  Half  of  the 
samples  received  a  pre-oxidation  anneal  in  nitrogen  at  the  oxidation 
temperature.  The  plot  of  log  oxide  thickness  versus  log  oxidation  time 
is  shown  in  Figure  3.  For  undoped  polycrystalline  silicon  the  differences 
observed  in  the  oxidation  rates  between  unannealed  and  pre-oxidation  annealed 
samples  are  mostly  within  the  accuracy  of  the  oxide  thickness  measurements. 
This  is  not  unexpected  because  x-ray  diffraction  analysis  of  the  undoped 
polycrystalline  silicon  film  as  deposited  by  LPCVD  at  615°C  shows  a  predomi¬ 
nant  grain  orientation  of  (110)  and  little  changes  are  observed  in  the 
structure  of  polycrystalline  silicon  after  annealing  at  temperatures  up  to 
1000°C  [1].  Analysis  of  the  oxide  growth  data  indicates  that  the  oxidation 
kinetics  of  the  film  resembles  that  of  (110)  oriented  single  crystal  siicon. 
The  parabolic  rate  constant  B  was  found  to  be  slightly  lower,  by  approximately 
25%,  than  that  of  the  (111)  single  crystal  silicon;  while  the  linear  rate 
constant  B/A  was  quite  similar  to  that  of  the  (111)  single  crystal  silicon. 
These  results  are  in  agreement  with  available  data  comparing  (111)  and  (110) 
lightly  doped  single  crystal  silicon  oxidation  kinetics  [2,3].  The  activation 
energies  for  B  (38  kcal/mole)  and  B/A  (35  kcal/mole)  are  the  same,  within  the 
experimental  accuracy,  as  for  the  B  and  B/A  respectively,  of  the  lightly  doped 
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(Ill)  single  crystal  silicon  in  this  lower  temperature  range.  The  graphs 
of  log  B  and  log  B/A  versus  reciprocal  tenpeature  are  shown  in  Figures  4  and 
5,  respectively. 

The  oxidation  of  doped  and  undoped  polycrystall ine  silicon  films  in  1  atm 
dry  02  ambient  at  800°,  900°,  and  1000°C  have  also  been  completed.  Annealed 
and  unannealed  samples  were  included.  All  wafers  were  pulled  in  argon  after 
each  oxidation  run.  This  completes  the  runs  of  the  atmospheric  part  of  the 
program.  Analysis  of  the  oxidation  kinetics  of  these  samples  and  conparison 
with  respect  to  single  crystal  silicon  oxidation  in  dry  02  are  underway. 

3.3  Structural  Evaluation 

Several  samples  of  doped  polycrystalline  silicon  oxidized  in  pyrogenic 
steam  at  1  atm  and  10  atm  have  been  submitted  for  cross-section  TEM 
evaluation.  Samples  were  chosen  to  evaluate  the  effects  of  pressure,  tem¬ 
perature,  and  annealing  on  oxides  of  the  same  thickness  range  (2000-2500  A) 
as  well  as  various  thicknesses  under  constant  temperature  and  pressure 
conditions. 
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PHOSPHORUS  DOPED  POLY-SI 
PYROGENIC  STERN  ,  1  RTM 
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FIGURE  1:  Residual  poly-silicon  film  resistivity  versus  poly-silicon 
oxide  thickness  grown  in  1  atm  pyrogenic  steam  at  750°, 
800°,  and  850°  C.  Results  for  annealed  and  unannealed 
phosphorus  doped  poly-silicon  film  are  shown. 


RESISTIVITY 


POLY-OXIDE  THICKNESS,  x0  </um> 

FIGURE  2:  Residual  poly-silicon  film  resistivity  versus  poly-silicon 
oxide  thickness  grown  in  10  atm  pyrogenic  steam  at  750°, 
800  ,  and  850°  C.  Results  for  annealed  and  unannealed 
phosphorus  doped  poly-silicon  film  are  shown. 
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FIGURE  5:  Linear  rate  constant  B/A  versus  1000/T  for  annealed  i 

and  unannealed  undoped  poly-silicon  films  oxidized  in 
pyrogenic  H?0  at  1  atm  and  750°  -  850°  C.  Linear  rate  j 

constants  of  (111)  Si  are  shown  for  comparison.  j 
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4.  OXIDATION  ENHANCEMENT  STUDIES 

W.  A.  TILLER,  Y.  T.  THATHACHAR I ,  W.  E.  DIBBLE,  JR., 
D.  N.  MODLIN  and  E.  M.  YOUNG 


4.1  Field  Enhancement 

In  the  field-enhanced  oxidation  studies,  the  native  oxide  thickness,  Xfj 
is  much  less  than  the  total  oxide  thickness,  Xg,  at  t=700°  C  with  a  <100> 
substrate  and  I  -  -  5pA.  Thus,  we  decided  to  check  to  see  if  the  experimental 
data  satisfied  a  constant  current-time  product  since  this  is  what  is  expected 
theoretically;  i.e., 

XB  XN  <No  (  } 

where  Xg  is  the  average  enhanced  oxide  thickness.  The  question  is,  does  this 
also  happen  for  the  peak  oxide  thickness,  i.e.,  do  we  have 

Xg  -  X^  a  It  (2) 

In  Figure  1,  we  see  the  superposed  results  of  an  I  =  -  5[iA  and  t  =  2.0  hr 
experiment  with  one  at  I  =  -  lO^A  and  t  =  1.0  hr.  Although  the  profiles  are 
quite  similar,  the  -  10nA  curve  is  wider  at  its  base  and  has  a  lower  peak 
height.  Thus,  altnough  eq.  1  holds,  eq.  2  does  not  hold  exactly.  This 
situation  is  readily  explained.  In  the  point  to  plane  discharge,  increasing  I 
leads  to  enhanced  natural  repulsion  betwen  the  ions  so  that  the  ion  beam 
widens  because  the  axial  electric  field  is  not  increased  in  proportion  to  the 
increased  space  charge.  Thus,  the  current  density  at  the  center  of  the  beam 
must  decrease  because  the  integral  of  the  total  beam  intensity  is  equal  to  the 
beam  current. 
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The  oxygen  partial  pressure  dependence,  Pg^,  was  investigated  with  the 
negative  ion  beam  at  700°  C  using  an  admixture  of  oxygen  in  argon.  The  values 
of  Xg  for  t  -  1  hr  at  I  =  -  IOjaA  with  10%,  50%,  80%,  and  100%  02  in  argon  are 
given  in  Figure  2.  We  note  that,  as  Pg^  is  reduced  by  a  factor  of  10,  Xg 
decreases  by  only  -  18%.  The  explanation  for  this  result  is  thought  to  depend 
upon  the  electron  attachment  probability  to  0  and  02. 

As  Pn  decreases,  the  concentration  of  unattached  electrons  entering  the 
u2 

oxide  increases.  We  expect  that  the  concentration  of  free  electrons,  n,  will 
be  given  by 

T)  -  Kx  expf-e^kT)  +  K2  expf-f^/kT)  (3) 

where  K  is  a  constant  and  e  is  the  binding  energy  of  an  electron  to  an  0  or  02 
(electron  affinity).  In  eq.  3,  K  will  decrease  as  Pg^  increases  and  will 
depend  upon  the  relative  population  of  0  and  02.  If  the  concentration  of  0 
stays  relatively  constant  even  though  the  concentration  of  02  decreases,  since 
e0  >  eg2,  one  might  expect  Xg  to  be  only  slightly  dependent  upon  Pg^.  Because 
of  eq.  3,  at  higher  T  for  a  given  t  and  I,  we  should  expect  Xg  to  fall  off. 
This  is  observed  experimentally. 

In  these  studies,  a  series  of  baseline  control  experiments  were  conducted 
The  Si  wafer  rested  on  a  quartz  plate  at  a  distance  of  0.5  cm  below  the  plati¬ 
num  needle  and  no  external  potential  was  applied  to  the  system.  However, 
directly  beneath  the  needle,  the  oxide  growth  was  enhanced.  This  was  later 
found  to  be  due  to  (i)  a  small  remnant  positive  polarization  potential  in  the 
quartz  plate  and  sapphire  standoffs  of  the  apparatus  and  (ii)  to  a  lesser 
extent  the  presence  of  platinum  in  the  system. 
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In  Figure  3,  two  sets  of  control  experiments  at  900°  C  with  Si  <1 1 1>  sub- 


trates  with  zero  applied  potential  are  shown.  The  first  set,  which  were 
performed  near  the  beginning  of  the  experimental  work,  reveal  an  enhanced  peak 
oxidation  rate  compared  to  that  expected  in  a  standard  oxidation  furnace.  The 
second  set,  conducted  some  months  later  with  no  changes  in  the  experimental 
apparatus,  reveal  an  even  greater  enhanced  peak  oxidation  rate  compared  to  the 
first  set.  It  was  anticipated  that  these  unwanted  electrical  effects  were 
associated  with  (i)  electrical  polarization  in  the  apparatus  due  to  the  flow 
of  current  I  and  (ii)  work  function  differences  between  dissimilar  materials; 
i.e. ,  the  equivalent  voltage  for  the  system,  Vg  ,  is  given  by 


'WF 


(4) 


We  expect  that  Vpol  >>  VWF  and  that  the  difference  between  set  1  and  set  2  of 
Figure  3  is  due  to  changes  in  Vp0] .  The  measured  volume  conductivity  at  900°C 
is  °v  «I0-9ohm~lcm~l  while  the  DC  resistance  for  a  fully  polarized  quartz 
plate  in  the  apparatus  is  -1.4  x  107  ohms  which  leads  to  an  IR  drop  of  »  70 
volts  across  the  quartz  plate.  Thus,  if  internal  polarization  completely 
cancels  this  voltage  drop  during  a  run,  the  remnant  field  left  in  the  system 
system  will  be  <!  -70  volts. 

Figure  4  shows  three  oxide  thickness  profiles  for  Si  <iii>  substrates 
at  900°  C  with  t  *  2  hrs.  Experiment  #63  was  an  open-circuit  control  while 
#119  and  #131  had  +50V  and  +100V,  respectively,  applied  to  the  needle.  As  #63 
was  the  first  experiment  performed,  the  oxidation  enhancement  is  indicative  of 
the  remnant  polarization  present  in  the  quartz  plate  at  that  time.  Experiment 
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#119  was  conducted  shortly  thereafter  and,  because  the  +50V  is  applied  in 
series  with  V  i  ,  the  resulting  thickness  profile  is  further  enhanced.  On  the 
other  hand,  experiment  #131  was  performed  after  a  nunber  of  experiments 
utilizing  a  beam  current  of  I  =  +5pA.  We  see  that,  for  an  applied  voltage  of 
+100V,  substantially  less  enhancement  is  obtained  than  that  found  in  experi¬ 
ment  #119.  This  is  because  Vp^changed  both  magnitude  and  sign  because  of  the 
positive  ion  current  forced  through  the  sample  and  the  quartz  plate. 

Assuming  a  linear  relationship  between  the  percentage  peak  enhancement 
and  the  resultant  voltage  on  the  needle,  one  can  estimate  the  voltage  due  to 
remnant  polarization.  Experiments  #63  and  #119  show  37%  and  62%  enhancement, 
respectively,  which  would  lead  to  a  V  -j  =  +74V  which  agrees  quite  well  with 
the  estimate  of  70  volts  obtained  earlier. 

4.2  Photo  Enhancement 

The  purpose  of  the  present  work  was  to  examine  the  extent  to  which  a 
true  photonic  or  optically  induced  enhancement  mechanism  is  at  work  in  silicon 
versus  enhancement  caused  by  local  excess  heating  of  the  surface  at  the  laser 
beam  spot.  In  our  study,  a  Lexel  argon-ion  laser  was  used  to  generate  a 
photon  flux  of  selectable  energy  and  adjustable  intensity.  In  the  continuous 
beam,  only  one  wavelength  of  the  five  most  powerful  lasing  wavelengths  was 
used  at  a  time;  i.e.,  2.4  eV,  2.5  eV,  2.54  eV,  2.6  eV  or  2.7  eV.  The  beam 
was  directed  into  a  thermal  oxidation  furnace  so  as  to  hit  the  center  surface 
of  a  vertical  sample  of  silicon  (see  Figure  5)  at  a  few  degrees  off  of  the 
normal  to  avoid  multiple-reflection  measurement  errors  of  the  beam  spot  power 
density.  The  silicon  sample  was  otherwise  exposed  to  typical  thermal  oxida¬ 
tion  processing. 
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Typical  oxidation  conditions  were  two  hours  at  900°  C  in  dry  02  which 
resulted  in  control  oxides  of  •  450  A  and  •  300  A  for  the  <11 1>  and  <100> 
orientations,  respectively.  The  peak  oxidation  enhancement,  measured  via 
ell ipsometry,  varied  linearly  with  beam  power  density  and  was  in  the  3%  to  30% 
range,  as  ilustrated  in  Figure  6.  We  note  that  the  enhancement  is  greater  for 
the  <100>  than  for  the  <1 1 1>  Si  which  is  opposite  to  that  expected  for  a 
purely  thermal  result. 

Experiments  were  also  carried  out  at  low  but  constant  beam  power  density 
while  the  photon  flux  was  varied  with  the  photon  energy  in  the  range  2.4  eV  to 
2.7  eV.  At  constant  beam  power  PQ  the  photon  flux  F  is  directly  proportional 
to  the  optical  wavelength  \  according  to 


where  c  is  the  speed  of  light  and  h  is  Planck's  constant.  Thus,  at  long  wave¬ 
lengths,  more  photons  arrive  at,  and  interact  with,  the  active  silicon  surface 
than  do  so  at  shorter  wavelengths.  Figure  7  illustrates  typical  results  for 
both  <1 1 1>  and  <100>  silicon  and  shows  the  percent  oxidation  enhancement  to 
increase  linearly  with  beam  wavelength  and  thus  with  photon  flux. 

A  variation  of  these  two  experiments  was  performed  to  conclusively 
demonstrate  the  nature  of  the  photonic  mechanism.  Here,  the  power  of  the 
longer  wavelength  beam  was  lowered  to  reduce  F^  to  the  same  value  as  that 

from  a  shorter  wavelength  beam  operating  at  Ps;  i.e. ,  we  require  that 
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If  a  photonic  mechanism  of  oxidation  enhancement  is  indeed  operative,  the 
amount  of  enhancement  produced  by  both  wavelengths  should  now  be  identical. 
The  results  shown  in  Figure  8  confirm  this  hypothesis. 

A  great  deal  of  care  must  be  exercised  over  the  course  of  a  typical  2 
hour  oxidation  experiment  both  in  measuring,  and  in  maintaining  constant,  the 
power  density  of  the  laser  beam.  This  is  a  prerequisite  to  observing  the 
photonic  effect  with  the  precision  of  the  experiments  becoming  more  critical 
as  the  wavelength  difference,  AX  *  |  Vj  -  Xj  | ,  decreases.  In  addition,  as 
the  beam  power  density  increases,  more  heating  of  the  silicon  surface  is 
likely  to  occur,  eventually  washing  out  the  photonic  effect.  Preliminary 
experiments  show  that,  as  the  beam  power  density  increases  for  fixed  ax, 
this  is  indeed  the  case. 
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thickness  vs  position  on  wafer  for  700  C  oxidations 


DISTRNCE  (INCHES) 


FIGURE  3:  Relationship  of  enhanced  oxidation  kinetic  results  to  "normal"  oxidation  data. 


900  C  < 1 1 1 >  2  Hrs 


FIGURE  4:  Oxide  thickness  profiles  for  900°  C,  <111>,  2  hr  oxidations. 
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FIGURE  5:  Experimental  set-up  for  laser  enhanced  oxidation  studies. 
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FIGURE  8:  Oxidation  enhancement  verses  laser  beam  power 
density  for  two  laser  wavelengths. 
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5.  Si /Si 0 2  INTERFACE  CHARGES 
A.  AKINWANDE,  C.  HO,  J.  PLUMMER 


In  the  last  annual  report,  it  was  demonstrated:  (a)  as  grown  charges 
follow  the  Deal  Np  triangle  qualitatively  but  a  careful  examination  of  the 
results  indicate  a  breakpoint  at  a  temp  of  about  950°  C.  This  led  us  to 
suggest  that  there  may  be  a  viscoelastic  relaxation  process  occurring  at  the 
interface  at  the  same  time  as  charge  generation  during  the  oxidation  process 
(Figure  1);  (b)  annealed  charges  appear  not  to  show  the  temperature  indepen¬ 
dence  that  the  bottom  part  of  the  Np  triangle  would  suggest.  There  seems  to 
be  a  break  point  at  about  950°  C  also.  This  led  us  to  suggest  that  the  anneal 
process  may  be  a  viscoelastic  decay  process  (Figure  2);  (c)  As  confirmation  of 
the  fact  that  the  annealed  charges  are  not  independent  of  temperature,  thermal 
cycling  experiments  were  carried  out  on  (111)  Si  wafers  oxidized  in  dry  02  at 
1000°  C  for  75  minutes.  These  wafers  were  later  annealed  in  argon  (<  2  ppm 
H20)  and  cycled  between  900°  C  and  1100°  C.  The  results  did  show  that  oxide 
charges  do  thermally  cycle  under  the  conditions  given  (Figure  3). 

Since  the  above  summarized  experiments,  numerous  experiments  have  been 
carried  out  to  answer  some  of  the  questions  raised  and  attempts  have  also  been 
made  to  develop  theoretical  tools  to  model  the  processes.  The  major 
accomplishments  over  the  last  few  months  are  itemized  below: 

(1)  A  fast  pull  technique  was  developed  to  enable  a  quench  of  Si 
wafers  from  processing  temperature  to  room  temperature.  This 
set  up  is  shown  in  Figure  4.  Quench  durations  of  <  3  secs  have 
generally  been  achieved. 
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(2)  In  the  first  set  of  experiments  using  the  new  experimental  set-up, in 
oxides  were  grown  at  1000°  C  in  dry  02  for  75  mins  and  were  fast 
pulled  (-  2  sec  pull).  They  were  subsequently  annealed  at  various 
tenperatures  (800°  C,  900°  C  and  1000°  C)  in  argon  which  has  <  1  ppm 
H20.  The  results  are  as  shown  in  Figures  5,  6  and  7.  From  these 
results,  the  anneal  process  appears  to  be  exponential  with  a  time 
constant  which  decreases  with  temperature.  The  anneal  process 
saturates  to  a  final  value  which  also  decreases  with  temperature. 

Finally,  this  has  led  to  a  suggestion  that  the  anneal  kinetics  may 
be  described  by: 

Nf(t)  =  (Nf(o)  -  Nfe(T))  exp  fV^T))  +  Nfg(T)  (1) 
where 

t(T)  =  K  exp(Ea/kT)  (2) 

In  order  to  explain  this  behavior,  we  favor  a  model  which  invokes  a 
maxwellian  viscoelastic  relaxation  process  at  the  Si /Si 0 2  interface. 

This  process  relieves  the  strain/stress  at  the  interface.  This 
implies  that  we  must  make  an  a  priori  assumption  that  strain/stress 
generated  by  the  lack  of  free  volume  during  the  oxidation  process  is 
the  origin  of  oxide  charges.  It  must  be  reiterated  though  that  this 
is  just  one  of  many  possible  models  and  careful  further  experimen¬ 
tation  is  required. 

We  have  also  investigated  a  situation  in  which  the  viscoelastic 
process  is  not  maxwellian.  This  would  result  in  a  situation  where 
the  fictive  "growth"  temperature  of  the  oxide  is  different  from  the 
anneal  temperature  [1].  It  is  expected  under  such  circumstances 
that  the  anneal  time  constant  would  be  a  function  of  the  processing 
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history  of  the  Si/$i02  wafers.  Preliminary  experiments  are  being 
carried  out  to  test  this  effect. 

(3)  The  second  set  of  recent  experiments  were  to  reproduce  the  thermal 
cycling  of  oxide  charges  under  different  anneal  conditions  (i.e. , 
temperature).  Oxides  were  grown  on  (111)  and  (100)  Si  Cz  wafers  at 
1000°  C  in  dry  02  for  75  mins.  They  were  subjected  to  20  min  pre 
heat  in  Argon  prior  to  oxidation,  20  min  1%  HCl/Ar  etch  +  20  min 
purge  in  Argon  at  1000°  C.  The  wafers  were  fast  pulled  from  the 
oxidizing  ambient  («  2  sec)  and  were  subsequently  given  anneals  in 
argon  and  cycled  between  800°  C  and  1100°  C.  As  the  results  shown 
in  Figure  8  indicate,  the  charges  do  thermally  cycle  in  argon  on 
both  orientations.  It  should  be  noted  that  the  annealed  charges  at 
800°  C  for  (111)  wafers  are  greater  than  the  as  grown  values  at 
1000°  C.  This  observation  coupled  with  the  fact  that  the  charges 
reproduce  under  thermal  cycling  has  led  to  the  suggestion  that  there 
may  be  a  reversible  reaction  going  on  at  the  interface  between  Si 
and  Si02,  that  results  in  an  equilibrium  value  of  oxide  charges 

for  a  particular  tenperature.  Attenpts  have  been  made  to  formulate 
this  problem  mathematically  and  results  thus  far  indicate  that 
more  careful  studies  are  required. 

(4)  The  third  set  of  new  experiments  were  to  study  the  as-grown  charge 
kinetics.  As  we  suggested  earlier,  it  has  been  suggested  in  the 
literature  [2]  there  may  be  a  possibility  that  during  an  oxidation 
process,  there  is  simultaneous  generation  and  annihilation  (anneal) 
of  oxide  charges  at  the  Si/Si02  interface.  The  problem  was  formu¬ 
lated  mathematically  under  the  assumptions  that:  (a)  there  is  a 
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distribution  of  oxide  charges  in  the  oxide  with  the  highest  value 
at  the  interface;  (b)  charges  are  generated  as  the  interface 
moves  past  the  point  in  question  in  the  oxide;  (c)  charges  are 
rather  immobile;  (d)  charges  at  a  point  in  the  oxide  anneal  with 
an  exponential  decay  starting  from  the  time  of  generation,  i.e., 
the  time  the  interface  moved  past  the  point  in  question.  The  time 
constant  would  be  what  is  expected  from  a  viscoelastic  relaxation 


t 


nil) 

=  G 


(3) 


and,  (e)  a  weighted  sum  of  the  charges  is  taken.  The  dependence 
this  model  would  predict  on  oxide  thickness  is  as  shown  in  Figure 
9.  Prel iminary  experiments  have  been  carried  out  and  the  results 
for  (111)  Si  wafers  are  shown  in  Figure  10.  It  is  our  belief  that 
by  extension  of  the  domain  of  the  experiment,  either  by  growing 
thicker  oxides  or  growing  oxides  at  higher  temperature  where  T  would 
be  shorter,  the  peak  predicted  by  the  model  would  be  observed. 
Experimentation  is  in  progress. 

(5)  The  fourth  set  of  experiments  were  to  investigate  the  as  grown 
charges  as  a  function  of  oxidation  temperature  at  constant  oxide 
thickness  (600 A) .  As  shown  in  Figure  11,  we  reproduce  earlier 
results  in  which  we  observe  a  breakpoint  at  a  temperature  of  about 
950°C. 


In  order  to  conclude  this  report,  it  can  be  seen  that  many  more  questions  have 
been  raised  than  have  been  answered.  We  have  thus  set  the  following  priori¬ 
ties: 
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(a)  Repeat  the  thermal  cycling  experiments  using  both  CZ  and  FZ  Si 
wafers.  This  would  eliminate  any  questions  as  to  the  effect  of 
dissolved  oxygen  on  Si  on  the  Si/Si02  interface. 

(b)  Study  the  anneal  kinetics  from  the  point  of  view  of  a 
reversible  process. 

(c)  Complete  the  study  of  the  effect  of  thermal  history  (i.e. , 
oxidation)  on  the  anneal  process. 

(d)  Determine  if  the  anneal  process  which  seems  to  be  limited  by  a 
bond  breaking/thermal  mechanism  can  be  aided  by  an  external  source 
of  energy  such  as  UV  light  at  low  terqperature  to  shorten  the  decay 
time  constant. 

(e)  Complete  the  experiments  on  the  generation  mechanisms. 
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FIGURE  1:  As-grown  (fast  pulled)  fixed  oxide  charge 
density  Nf  vs.  growth  temperature  for  (ill) 
and  (100)  wafers.  While  the  general  behavior 
expected  from  the  Nf  triangle  is  observed,  a 
distinct  break  also  occurs  at  ^950°  C. 


FIXED  OXIDE  CHARGE  Np  (cnf2) 


FIGURE  3:  Fixed  oxide  charge  density  Nf  vs.  Argon  anneal 
cycle  for  (111)  wafers  thermally  cycled  between 
900°  C  and  1100°  C  in  Ar,  This  illustrates  that 
the  bottom  part  of  the  Nf  triangle  is  not  flat. 
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FIGURE  4:  Schematic  representation  of  fast  pull  apparatus  designed  to  permit  <3  sec 
quench  times  from  oxidation  temperature  to  room  temperature. 
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FIGURE  5: 


Example  of  anneal  kinetics  of  oxide  charges  for  (100) 
wafers  in  Argon  at  800°  C  and  1000°  C.  Note  that  both 
decay  time  constants  and  the  final  annealed  values  are 
different  at  the  two  temperatures. 
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OXIDE  THICKNESS  Cflna) 

FIGURE  9:  Predicted  dependence  of  on  oxide  thickness  based 
upon  the  model  described  in  the  text,  and  a  visco¬ 
elastic  relaxation  time  of  4  hours. 
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Experimental  dependence  of  Nf  oxide  FIGURE  11:  As  grown  Nf  charge  densities  as  a  function 
thickness  for  wafers  oxidized  at  1000°  C  of  oxidatiin  T^for  oxide  grown  to  a  simila 
in  dry  0?.  thickness  (600^)  at  each  T.  Note  the 


6.  MEASUREMENT  AND  CHARACTERIZATION  OF  THIN  OXIDES  ON  POLYSILICON 


E.  KUTLU,  K.  SARASWAT,  J.  PLUMMER 


Our  first  experiments  on  the  oxidation  of  plysilicon  showed  that  the 
oxidation  kinetics  of  poly  are  different  from  the  oxidation  kinetics  of  single 
crystal  silicon  mostly  during  the  linear  regime  [1].  It  is,  therefore, 
necessary  to  work  with  thin  oxides  in  order  to  understand  the  oxidation 
kinetics,  of  polysilicon.  Although  several  measurement  techniques  are 
available  for  the  measurement  of  thin  oxides  on  single  crystal  silicon,  they 
are  not  generally  useful  for  the  measurement  of  thin  oxides  on  poly. 

The  sensitivity  of  alpha  step  measurements  to  variations  of  the  oxide 
thickness  on  the  wafer  makes  the  accuracy  of  the  measurements  questionable. 

The  measurements  give  ±  20%  thickness  variations  for  oxides  thinner  than 
1000  A.  Eppl ipsometry  measurements  would  not  be  meaningful  because  of  the 
multiple  reflections  coming  from  the  multiple  layer  structure.  The 
spectrophotometer  technique  would  also  fail,  since  short  wavelengths,  which 
are  necessary  for  the  measurement  of  thin  films,  would  be  absorbed  by  the 
silicon. 

An  alternative  technique,  which  we  are  investigating,  uses  C-V 
measurements.  When  the  oxide  capaticance  is  known,  oxide  thickness  can  be 
obtained  from  the  relationship 
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where  t0x  is  the  oxide  thickness,  e0  is  the  permittivity  of  free  space, 

Kox  is  the  dielectric  constant  of  the  oxide  and  Cox  is  the  oxide  capacitance. 
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The  capacitor  structure  prepared  for  these  measurements  is  shown  in 
Figure  1.  In  the  following  sections  analyses  of  the  different  measurements  on 
this  structure  will  be  discussed,  and  a  set  of  experiments  which  are  in 
progress  will  be  described.  Although  the  structure  is  distributed  in  nature, 
the  equivalent  circuit  elements  will  be  treated  as  lumped  circuit  elements. 

6.1  Analysis  of  the  Measurements 

With  the  structure  shown  in  Figure  1,  three  different  measurements  seem 
possible:  substrate-ring  shorted;  substrate  floating,  and  a  2-step 
measurement  which  will  be  described  later. 

1.  Substrate-Ring  Shorted 

The  structure,  its  lump  equivalent  and  measured  equivalent  circuits 
are  shown  in  Figures  2a,  2b,  2c  respectively.  In  these  figures:  Cox 
is  the  capacitance  of  the  oxide  layer  on  the  silicon,  Cs  is  the  deple¬ 
tion  or  inversion  layer  capacitance  of  the  substrate,  RSp  the  spreading 
resistance  of  the  poly.  Rs  is  the  substrate  resistance  of  the  silicon 
and  Cm  and  gm  are  the  measured  equivalent  capacitance  and  conductance 
of  the  system  repectively. 

As  can  be  seen  from  Figure  2,  the  measured  capacitance  is  in  general 
different  from  the  oxide  capacitance  Cox,  depending  on  the  values  of 
the  parasitic  elements.  If  some  of  those  such  as  Rs,  Cs  and  C0  are 
known,  Cox  and  Rsp  can  be  found  from  Cra  and  gm  uniquely  for  the 
given  structure.  The  resultant  expressions  are  rather  complicated, 
however,  the  dependence  of  Cm/Cox  on  Rsp»  ^o  and  R  which  is  the 
radius  of  the  dot  are  given  in  Figures  3-6.  In  these  figures  it  is 
assumed  that  Rs  *  50Q  and  Cs  is  negligible. 
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As  can  be  seen  from  Figure  3  and  Figure  2b,  there  are  two  limit¬ 
ing  cases  corresponding  to  highly  and  lightly  doped  poly: 

a)  lim  r  =  r  (2) 

R  +  o  ^m  “  ^ox 

sp 


which  means  that  for  highly  doped  samples,  the  substrate  parasitic 
elements  t0,  Rs,  Cs  do  not  have  an  effect  on  the  measurements. 


b) 


1  im.  C 


R  ♦ 
sp 


ox 

2,  2r,  2 


1  +  “  Cox^Rs 


(3) 


where  =  (4) 

OX  ^0  OX  0  0  s 

which  means  that  depending  on  the  oxide  thickness  on  the  substrate,  the 

measured  capacitance  will  differ  from  the  true  oxide  capacitance. 

Neglecting  the  effect  of  the  depletion  capacitance  i.e.  Cq  »  C0, 

three  special  thicknesses  are  distinguishable. 


In  this  case,  under  the  condition  w2c£x2Rs2  <<  1  Cm  = 

For  two  different  thicknesses,  this  case  is  shown  in  Figure  5  and 
Figure  6. 

b2'  co  »  lox 

which  means  C0  «  Cox.  In  this  case,  since  C0  will  dominate,  we 
will  measure  C0,  rather  than  Cox. 
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to<< 


t 


ox 


or  Cox  <<  C0.  Since  we  will  measure  the  oxide  capacitance  directly, 
this  is  the  condition  we  have  to  consider  when  we  want  to  measure  thin 
oxides  on  lightly  doped  samples.  The  effect  of  t0  on  the  measured 
capacitance  for  lightly  doped  poly  is  shown  in  Figures  7-8. 

Now,  the  question  is,  in  the  case  of  lightly  doped  poly,  how  thin  can 
the  oxide  layer  on  the  substrate  be  made  in  order  to  get  the  correct  oxide 
thickness  on  poly?  We  cannot  avoid  having  an  oxide  layer  on  the  Si  sub¬ 
strate,  since  this  would  totally  change  the  polysilicon  oxidation  kinetics. 
Furthermore,  if  we  make  t0  thinner  than  200A,  then  its  thickness  becomes 
comparable  with  that  of  the  depletion  layer  and  we  can  no  longer  neglect 
the  depletion  capacitance.  Therefore,  with  this  measurement  technique, 
especially  for  the  lightly  and  moderately  doped  poly,  we  cannot  avoid  the 
effects  of  parasitic  elements  on  the  measurement. 

2.  Floating  Substrate 

Another  possible  measurement  which  can  be  done  between  the  ring  and 
dot  by  letting  the  substrate  float,  seems  to  avoid  the  parasitic  elements 
of  the  structure.  The  structure,  and  as  a  first  approximation  its  equiv¬ 
alent  circuit  are  shown  in  Figures  9a-9b  respecti vely.  Figure  9c 
illustrates  the  measured  equivalent.  As  is  seen  from  the  figures  the 
calculation  of  the  oxide  capacitance  Cox,  requires  a  simple  unique 
conversion  from  the  measured  capacitance  Cm  and  the  conductance  gm. 


ox 


■  (1  *  4-) 


Rsp’ 


U1 


K?) 


where 


2-2 
“  Cm 


(5) 
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Although  this  technique  looks  simple,  in  the  case  of  leakage  from 
the  poly  to  the  substrate,  the  problems  of  the  first  technique  would 
still  exist  here.  Furthermore,  since  the  whole  structure  is  distributed, 
a  first  order  lunped  circuit  approximation  might  not  be  adequate. 

3.  2-Step  Measurement 

This  technique  utilizes  a  2-step  measurement,  one  between  the  gate 
and  the  substrate,  and  the  other  between  the  ring  and  the  substrate. 

The  advantage  of  this  measurement  is  that  it  provides  knowledge  of  the 
parasitics  of  the  system.  The  second  measurement  will  give  directly  the 
parasitics  of  the  system,  and  by  using  this  information  in  the  first 
measurement,  we  will  be  able  to  extract  the  oxide  capacitance.  The 
measurement  structures  and  their  equivalent  circuits  are  shown  in  Figures 
10-11.  In  these  figures:  Cx  is  the  total  parasitic  capacitance  of  the 
system,  Rx  is  the  total  parasitic  resistance  of  the  system,  Rc  is  the 
contact  resistance  between  the  A1  and  the  poly,  and  Cox  is  the  oxide 
capacitance.  Cmi,  and  Cm2  are  the  measured  equivalent  capaci¬ 
tances  and  conductances  of  the  2-step  measurement.  From  Figures  10-11, 
it  can  be  shown  that  the  above  parameters  can  be  calculated  in  terms  of 
the  measured  capacitances  and  conductances  as  follows: 


Rc  = 


W1 


+  4) 


Rx  = 


W1 


' — T  '  Rc 


Cx  •  W1  ♦  3  ’ 


c  ■  [?rr  (: 


'm2  1  +  K; 


r)  -  f]  (6) 
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(7) 


where: 


2-2 
w  Cml 


dml 


“24 
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Experiments 

In  order  to  verify  the  above  analysis,  the  experiments  should  reflect  the 
effects  of  the  parasitic  elements  of  the  structure.  So,  the  experiments 
should  be  done  on  both  highly  and  lightly  doped  samples  which  have  thin  and 
thick  oxide  layers  on  the  substrte.  The  substrate  should  be  highly  doped  in 
order  to  reduce  the  ffect  of  Rs  and  Cs.  The  measurements  should  also  be 
done  for  different  thicknesses  of  oxide  grown  on  the  poly.  Given  the  above 
considerations,  the  following  samples  were  prepared. 

1.  Four  2",  <111>,  B  doped  wafers  with  0.01  Qcm  resistivity  were  used  as 
substrates. 

2.  Two  of  these  wafers  were  oxidized  in  dry  02  at  850°C  for  15  min  giving  an 
average  oxide  thickness  of  250A.  The  other  two  wafers  were  oxidized  in  dry  02 
at  1100°C  for  70  min  giving  an  average  oxide  thickness  of  1450A. 

3.  Poly  was  deposited  on  all  of  the  wafers  with  LPCVD.  The  average  thickness 
was  0.5(im. 

4.  The  halves  of  the  first  and  second  group  of  wafers  were  implanted  with  P 
with  a  dose  of  1015cm"2,  whereas  the  other  halves  were  implanted  with  a  dose  of 
1016cm-2. 

5.  All  of  the  samples  were  annealed  in  Ar  at  1000°C  for  70  min,  in  order  to 
allow  the  dopant  to  redistribute. 

6.  One-half  wafer  from  each  group  with  1015  doping  was  oxidized  at  850°C  for 
15  min.  Both  of  these  oxidation  conditions  should  give  200A  average  oxide 
thickness  on  the  poly.  The  rest  of  the  samples  were  oxidized  in  dry  02  at 
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1100°C  for  70  min.  This  last  oxidation  should  give  1500A  average  oxide 
thickness  on  the  poly.  Measurements  are  underway  at  the  present  time  to 
characterize  these  samples.  When  all  the  measurements  are  completed,  we  would 
expect  to  see  the  following  results. 

1.  With  the  floating  substrate  and  the  2-step  mesurements  we  would 
expect  to  measure  the  capacitance  of  the  oxide  on  all  of  the 
samples,  regardless  of  their  preparation. 

2.  With  the  shorted  substrate  measurement,  for  highly  doped  samples 

we  would  expect  to  measure  the  correct  oxide  capacitance.  However, 
for  lightly  doped  samples,  depending  on  the  oxide  thicknesses  both 
on  the  substrate  and  on  poly  we  should  observe  differences  on  the 
measured  values. 

The  development  of  this  technique  for  measuring  thin  oxides  on  polysilicon 
should  allow  us  to  characterize  the  growth  kinetics  of  SiO 2  on  this  film  (and 
also  on  silicides  and  other  thin  films)  and  hence  should  provide  better  data 
and  models  for  incorporation  in  SUPREM  III. 
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7.  ION  IMPLANTATION  AND  RANGE  STATISTICS 
M.  GILES,  L.  CHRISTEL,  J.  GIBBONS 

During  this  period  the  main  enphasis  of  our  work  has  been  the  application 
of  the  Boltzmann  transport  equation  to  the  calculation  of  ion  implantation 
profiles  for  light  ions  and  at  low  energies.  As  devices  shrink  to  sub-micron 
dimensions  the  thickness  of  overlying  films,  such  as  oxides,  must  be  corres¬ 
pondingly  reduced.  This  will  in  time  require  implantation  energies  to  be 
lowered.  It  is,  therefore,  necessary  to  develop  methods  of  calculating 
implanted  ion  profiles  at  these  low  energies.  Some  time  has  also  been  spent 
in  developing  a  very  simple  model  of  the  diffusion  of  implanted  ions  based  on 
a  joined  half-gaussian  profile,  suitable  for  hand  calculation  or  as  a  rapid 
approximation  in  the  early  stages  of  process  development. 

7.1  Low  Energy  Light  Ion  Implantation 

The  Boltzmann  transport  equation  approach  [1]  has  been  successfully 
applied  to  many  ion  implantation  problems,  predicting  range  and  damage 
distributions,  recoil  distributions  [2]  and  stoichiometric  disturbances  in 
compound  semiconductors  [3].  In  this  approach,  numerical  integration  of  the 
transport  equation  proceeds  stepwise  into  the  target.  At  each  step  the  state 
is  described  by  a  matrix  in  which  the  value  of  an  element  corresponds  to  the 
number  of  ions  with  a  certain  energy  E;  moving  at  an  angle  ©j  to  the  surface 
normal.  During  each  step  scattering  redistributes  ions  in  the  energy-angle 
matrix.  Ions  which  are  scattered  to  angles  greater  than  90°  are  called 
backscattered  and  were  previously  considered  stopped  at  the  scattering  point. 

However,  when  the  incident  Ion  is  much  lighter  than  the  target  atoms,  a  large 
fraction  of  the  incident  ions  are  backscattered.  To  follow  the  motion  of 
these  ions  a  multiple  pass  approach  has  been  developed,  where  each  pass 

follows  the  motion  of  ions  backscattered  during  the  previous  pass.  t 
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In  the  present  calculations,  backscattered  ions  are  removed  from  the 
matrix  and  their  energy  and  angular  distribution  is  stored  for  future  use.  At 
the  end  of  the  first  pass,  we  then  have  a  concentration  profile  as  before, 
together  with  a  set  of  matrices  of  backscattered  ions  for  each  step  into  the 
target.  The  second  pass  is  made  from  inside  the  target  towards  the  surface, 
which  is  the  direction  of  motion  of  the  previously  backscattered  ions.  At  the 
beginning  of  each  step  the  stored  backscatter  distribution  for  the  depth  is 
added  to  the  continuing  distribution  of  ions  still  moving  from  the  previous 
step.  Ions  backscattered  during  the  step  are  again  stored  in  preparation  for 
a  third  pass,  and  so  on  until  all  of  the  ions  have  stopped. 

The  calculated  profile  for  Boron  implanted  into  Silicon  at  15  keV  is 
shown  in  Figure  1,  where  the  target/ion  mass  ratio  is  2.55.  Here  we  observe 
that,  even  though  60%  of  the  incident  Boron  ions  are  backscattered  during  the 
first  pass,  there  is  little  change  in  the  final  profile.  This  is  because  most 
of  the  backscattered  ions  are  moving  almost  perpendicular  to  the  surface  so 
are  likely  to  be  stopped  or  backscattered  again  with  only  a  small  change. 
Similar  results  are  obtained  with  higher  beam  energies,  leading  to  the 
conclusion  that  single  pass  Boltzmann  or  LSS  [4]  techniques  will  give  satis¬ 
factory  profile  descriptions  for  these  cases. 

In  two  layer  targets,  different  effects  are  observed  depending  on  the 
relative  masses  of  the  target  atoms  In  the  two  layers.  When  the  layers  are  of 
similar  mass,  such  as  Silicon  Dioxide  on  Silicon,  there  is  no  great  change  in 
profile  except  for  some  Increase  In  recoil  mixing  at  the  Interface.  When  the 
layers  are  of  different  mass,  implanted  Ions  tend  to  recoil  from  the  heavier 
layer  and  accumulate  In  the  lighter  layer.  This  Is  Illustrated  by  Figure  2, 
which  shows  the  profiles  for  15  keV  Beryllium  Implanted  into  300A  Silicon 


Nitride  on  Gallium  Arsenide.  The  dose  retained  in  the  GaAs  is  reduced  from 
88%  to  75%  of  the  incident  beam,  a  very  substantial  reduction,  and  recoil 
mixing  is  increased. 

7 .2  Approximation  of  Implanted  Profiles  by  Joined  Half-Gaussians 

It  is  well  known  that  non- symmetrical  impurity  profiles  are  obtained  when 
common  ions  such  as  Boron,  Phosphorus  and  Arsenic  are  Implanted  into  amorphous 
or  randomly  oriented  Silicon.  It  has  been  shown  [5,6]  that  Pearson  distri¬ 
bution  functions  can  provide  an  accurate  fit  for  Impurity  distributions  for 
implantation  energies  up  to  several  hundred  keV.  Gibbons  and  Itylroie  [7] 
showed  that  a  distribution  consisting  of  two  half-Gaussian  distributions 
joined  at  a  modal  projected  range  could  provide  a  satisfactory  fit  to 
impurity  profiles  obtained  under  typical  implantation  conditions. 

We  have  chosen  to  construct  joined  half-Gaussian  distributions  which  have 
the  same  mode,  mean  and  standard  deviation  as  the  Pearson  distribution  for  the 
implanted  profile.  The  distributions  can  be  constructed  from  the  first  three 
moments  of  the  profile,  which  are  available  is  tabular  form  for  common  ion/ 
target  combination  [8], 

In  most  processes  ion  implantation  is  followed  by  several  high 
temperature  steps  during  which  the  implanted  distribution  diffuses.  The 
diffused  distribution  can  also  be  approximated  by  a  joined  half-Gaussian  with 
new  values  for  R^,  oi,  02*  Expressions  have  been  derived  for  these  para¬ 
meters  as  a  function  of  diffusivity  and  diffusion  time.  As  an  example.  Figure 
3  conpares  the  joined  half-Gaussian  approximations  with  a  numerical  solution 
for  the  diffusion  of  a  Pearson  distribution,  for  50  keV  Boron  in  Silicon.  It 
can  be  seen  that  the  accuracy  of  the  approximation  Improves  with  diffusion  as 
both  profiles  tend  towards  a  symmetric  Gaussian.  Considering  the  quality  of 
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fit  and  the  extreme  simplicity  with  which  Gaussian  distributions  can  be 
constructed,  it  seems  likely  that  the  joined  hal f-Gaussian  approximation  can 
provide  satisfactory  impurity  profile  estimations  for  first  order  analysis  of 
experimental  low  energy  implanted  distributions  and  provide  the  means  for 
calculating  profile  changes  produced  by  diffusion. 
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1:  implanted  profile  for  15  keV  Boron  into  Silicon  at  a  dose  10  cm' 
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CIGUR£  2:  Implanted  profile  for  15  keV  Beryllium  into  300  A  of  Silicon  Nitride  on 
Gallium  Arsenide  at  a  dose  1015cm2,  including  recoil  profiles  of  Silicon 
Nitrogen. 
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FIGURE  3:  Comparison  of  Pearson  IV  with  joined  half-Gaussian  distributions  for 
diffusion  of  Boron  implanted  at  50  keV  into  Silicon. 


8.  SILICIDE/SILICQN  CONTACTS 
F.  C.  SHONE,  K.  SARASWAT,  J.  PLUMMER 


During  the  last  few  years  ohmic  contacts  formed  early  in  the  process 
sequence  to  get  denser  device  structure  have  become  important.  Such  contacts 
are  subjected  to  several  high  temperature  steps  after  deposition.  This  modifies 
the  requirements  imposed  on  the  materials  being  used  to  form  the  contacts.  It 
has  been  observed  that  silicides  of  W,  Mo,  Ta,  Ti  are  compatible  with  IC  fabri¬ 
cation  technology.  They  have  fairly  high  conductivity,  they  can  withstand  all 
the  chemicals  encountered  during  fabrication  processes,  and  thermal  oxidation 
of  their  silicides  can  be  accomplished  in  oxygen  and  steam  to  produce  a  passi¬ 
vating  layer  of  Si 02 -  Their  contacts  to  shallow  p-n  junctions  are  reliable,  and 
fine  lines  can  be  etched  by  plasma  etching  these  materials.  This  indicates  that 
the  silicides  of  W,  Mo,  Ta,  and  Ti  will  be  used  as  interconnection  layers  and 
ohmic  contact  materials  in  future  VLSI  devices.  It  is  therefore  of  interest  to 
model  these  parameters  in  SUPREM  III. 

Mochizuki  et.al.  annealed  MoSi2  contacts  to  heavily  phosphorus  doped  N+ 
diffusions  and  observed  above  700°  C  the  contact  resistance  increased  rapidly. 
The  problem  was  caused  by  dopant  diffusion  from  the  silicon  to  the  silicide 
causing  depletion  of  donor  inpurities  at  the  silicon  surface  and  hence  resulting 
in  increased  contact  resistance.  To  overcome  this  problem,  they  doped  the 
silicide  heavily  with  phosphorus  prior  to  the  anneal  and  found  no  increase  in 
contact  resistance  took  place  even  after  annealing  up  to  1100°  C.  Obviously,  in 
this  case  the  loss  of  dopant  did  not  take  place  and,  therefore,  the  contact 
resistance  remained  stable. 

The  dopant  concentration  behavior  in  silicide  layers  depends  on  the  type  of 
dopant,  annealing  conditions,  which  films  are  present  and  the  stoichiometry  of 
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the  sil icicle.  The  diffusivity  of  the  dopants  has  been  found  to  be  more  than 
three  orders  of  magnitude  higher  than  in  single  crystal  silicon.  This  rapid 
diffusion  and  redistribution  can  explain  observed  changes  in  contact  resistance 
during  high  temperature  cycling.  If  these  sil icide  are  to  be  used  as  contact 
materials,  the  modeling  of  dopant  diffusion  in  silicires  and  segregation  and 
redistribution  at  the  sil icide/sil icon  interface  are  of  extreme  importance. 

In  experiments  conducted  in  our  laboratory,  we  have  deposited  WSi2  thin 
films  of  about  2000  A  thickness  on  silicon  substrates  and  annealed  then  in  N2 
for  30  minutes  at  900°  C  to  reduce  the  resistivity.  After  annealing,  we  heavily 
implanted  the  film  with  As+  ions  and  then  put  CVD  oxide  on  it  as  a  cap. 

Finally,  the  wafers  were  placed  in  a  furnace  to  diffuse  the  impurities  into  the 
silicon  substrate. 

We  used  RBS  to  measure  the  impurity  redistributions  in  the  sil icide  layer 
spreading  resistance  to  measure  the  inpurities  profile  in  the  silicon  and  AES  to 
identify  the  interface  segregation  behavior.  It  was  found  that  the  As+  dopant 
piled  up  at  the  interface  and  the  dopant  profile  in  silicon  was  similar  to  that 
found  with  doped  polysilicon  diffusion  sources. 

The  initial  experimental  results  show  that  the  doped  sil icide  acts  as  a 
diffusion  urce  to  form  N+  junctions  and  simultaneously  forms  an  ohmic  contact 
quite  easily.  Further  study  will  emphazise  the  dopant  segregation  behavior  and 
the  mechanisms  of  dopant  diffusion  in  the  sil icide  and  the  silicon  substrate  to 
develop  quantitative  models  suitable  for  SUPREM  III  implement. 


9.  GETTERING  AND  TRANSIENT  PROCESS  CHARACTERIZATION 
G.  BRONNER,  C.  HO,  J.  PLUMMER 

The  past  work  of  this  project  has  shown  that  tungsten  silicide  (WSi 2)  can 
be  used  as  a  metal  for  contacting  devices  without  degrading  device  performance. 
Wsi2  can  withstand  high  temperatures ,  one  can  anneal  a  finished  device  at  high 
temperature  (anywhere  from  800°  C  to  1200°  C),  cool  the  device  back  to  room 
temperature,  and  then  measure  how  the  anneal  has  affected  the  device  character¬ 
istics.  Thus  by  doing  a  series  of  anneals  and  device  mesurements,  a  transient 
process  can  be  characterized.  This  technique  has  been  applied  to  the  problem  of 
backside  damage  gettering.  Bipolar  transistors  were  built.  By  looking  at  their 
low  current  behavior  one  can  infer  the  level  of  metallic  contamination  in  the 
device.  Studying  gettering  with  this  technique  raised  a  number  of  questions 
that  we  have  been  looking  at  over  the  past  few  months. 

The  first  question  we  asked  was  whether  we  were  really  looking  at  metallic 
contamination  when  we  examined  the  behavior  of  the  base  current  in  bipolar 
transistors.  It  is  well  known  that  metals  can  affect  bipolar  device  character¬ 
istics;  the  issue  was  whether  they  were  the  only  things  that  could  be  affecting 
the  devices.  Recently,  it  has  been  suggested  that  vacancies  can  affect  device 
performance  [1].  This  seems  reasonable  on  several  levels.  Vacancies  are  known 
to  affect  many  processes  in  silicon.  In  SUPREM,  diffusion  and  oxidation  in 
heavily  doped  regions  are  explained  via  vacancy  arguments.  Additionally, 
vacancies  are  present  in  the  right  numbers.  At  850°  C  the  total  number  of 
vacancies  in  lightly  doped  silicon  is  8.6  x  10"  cnr3.  In  doped  regions  these 
values  are  even  higher.  If  these  vacancy  levels  could  be  frozen  in  and  if  they 
were  electrically  active,  they  could  swamp  out  any  effects  due  to  metallic 
contamination. 


-79- 


To  test  this  idea,  the  following  experiment  was  done.  A  wafer  was  oxidized 
at  1000°  C  to  passivate  the  surface.  It  was  then  annealed  at  temperatures 
between  600°  C  and  1000°  C.  After  each  anneal  the  wafer  was  cooled  as  quickly 
as  possible  to  room  temperature.  It  then  received  a  1/2  hour  anneal  at  450°  C 
in  forming  gas.  Bulk  lifetime  was  then  measured  by  looking  at  the  change  in 
conductivity  as  a  light  was  strobed  on  and  off. 

If  frozen-in  vacancies  limit  device  performance,  then  lifetime  should  drop 
off  with  higher  anneal  temperature.  A  plot  of  lifetime  vs.  inverse  temperature 
(t  vs.  1/T)  should  give  the  activation  energy  of  vacancy  formation.  The  data  is 
plotted  in  this  form  in  Figure  1.  We  see  that  from  700°  C  to  1000°  C  the 
life-time  is  essentially  constant.  Only  at  the  lowest  anneal  temperature, 

600°  C,  is  the  lifetime  somewhat  higher.  From  this  we  conclude  that  frozen-in 
vacancies  do  not  seem  to  directly  affect  device  behavior.  Thus,  when  we  look 
at  the  low  current  operation  of  bipolar  transistors  we  are  likely  seeing  a 
phenomena  caused  by  metallic  contamination. 

The  reason  for  the  rise  in  lifetime  with  the  600°  C  anneal  is  not  clear, 
but  it  is  possible  to  hypothesize  an  explanation.  When  metals  are  in  solution 
in  si  icon,  they  can  either  be  on  or  off  lattice  sites,  i.e. ,  they  can  be  substi¬ 
tutional  or  interstitial.  In  general,  interstitial  diffusion  is  much  faster 
than  substitutional  diffusion  for  these  metals.  Once  equilibrium  is  reached, 
solubilities  tend  to  be  much  higher  for  the  substitutional  forms  of  the  metal. 
Only  the  substitutional  metals  are  electrically  active  in  silicon.  Thus,  when 
metals  are  first  introduced  into  silicon,  they  are  interstitial.  These  inter¬ 
stitial  atoms  diffuse  quickly  until  equilibrium  is  established.  They  then  would 
Mke  to  drop  onto  substitutional  lattice  sites.  For  that  to  happen  the 
following  reaction  must  take  place: 

interstitial  +  V  institutional ,  (I) 
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where  M  stands  for  a  metal  atom  and  V  stands  for  a  silicon  vacancy.  This 
reaction  is  a  simplified  one,  but  it  does  describe  qualitatively  what  takes 
place.  It  has  recently  been  suggested  that  instead  of  a  substitutional  metal 
atom  one  might  have  a  vacancy-metal  interstitital  complex  [2,3],  This  effect, 
in  addition  to  the  various  charge  states  of  the  atoms,  is  ignored. 

For  this  reaction  we  can  write  the  affinity  product. 


J 


C^sub 
C^int  C 


(2) 


where  C^sub  is  the  concentration  of  substitutional  metal  atoms,  C^int  is  the 
concentration  of  interstitial  metal  atoms,  and  C y  is  the  concentration  of 
vacancies.  When  J  is  less  than  the  equilibrium  constant  k,  the  reaction  is 
driven  to  the  right,  i.e. ,  all  the  metal  atoms  assume  substitutional  sites. 

This  is  the  situation  at  high  temperatures  when  Cy  is  relatively  high.  As 
the  temperature  of  the  anneal  decreases,  Cy  decreases  and  J  increases.  When  J 
becomes  greater  than  k,  the  reaction  tends  to  reverse.  Substitutional  gold 
jumps  off  of  the  lattice  to  form  a  vacancy  and  an  interstitial  metal  atom. 
Evidently,  at  600°  C  the  concentration  of  vacancies  drops  enough  to  cause  this 
to  happen.  Thus,  a  fraction  of  the  metal  atoms  become  interstitial  (and  also 
electrically  inactive)  which  we  see  as  an  increase  in  silicon  lifetime. 

With  this  background  we  feel  that  our  experiments,  with  anneals  in  the 
range  of  800°  C  to  1000°  C,  do  indeed  look  at  the  metal  concentrations  in  our 
devices.  Therefore,  we  have  fabricated  a  number  of  wafers  with  various  test 
structures.  These  include  different  sized  bipolar  transistors,  M0S  capacitors, 
and  various  contact  resistance  structures.  The  bipolar  transistors  allow  us  to 
indirectly  follow  the  levels  of  metallic  contamination  by  looking  at  their  low 
current  operation.  Different  sized  transistors  let  us  separate  the  bulk 
effects,  which  we  are  interested  in,  from  any  surface  effects.  The  MOS 
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capacitors,  through  pulsed  lifetime  tests,  provide  an  independent  measure  of 
metallic  contamination.  The  contact  resistance  structures  look  at  the  Wsl2/Si 
interface.  From  these,  we  will  be  able  to  see  the  amount  of  dopant  segregation, 
and  determine  whether  any  changes  in  device  characteristics  are  due  to  these 
effects.  All  of  these  devices  have  been  built  in  float  zone  material  to 
eliminate  competition  between  the  backside  getters  being  studied  and  internal 
gettering  due  to  oxygen  precipitation. 

Before  following  the  time  behavior  of  the  backside  getter,  there  are  a  few 
facts  we  were  interested  in  establishing.  The  first  of  these  had  to  do  with  how 
WSi 2  contacts  affect  device  performance.  All  of  our  previous  work  has  shown 
that  silicide  contacts  do  not  degrade  performance.  There  is  the  possibility 
that  they  could  improve  performance.  A  recent  paper  by  Thompson  and  Tu  [4]  has 
shown  that  aluminum  contacts  on  silicon  can  getter  copper,  silver,  and  gold.  If 
WSi 2  acts  in  a  similar  fashion,  we  might  have  trouble  following  the  behavior  of 
a  backside  getter. 

To  see  how  much  of  a  problem  this  might  be,  we  have  designed  some  control 
experiments  to  test  the  limits  of  our  experimental  procedure.  The  first  of 
these  looks  at  the  behavior  of  a  backside  getter  with  and  without  WSi 2 
contacts.  A  set  of  wafers  consisting  of  the  test  structures  outlined  earlier  is 
fabricated  up  to  the  metal  deposition  step.  One  wafer  is  taken,  receives  a 
backside  getter,  and  then  an  8  hour  gettering  anneal.  Contacts  are  made  to 
devices  using  Aluminum  and  the  devices  are  characterized. 

Another  wafer  receives  a  backside  getter,  but  before  receiving  a  getter 
anneal,  its  devices  are  contacted  using  WSi 2*  It  then  is  annealed  for  8  hours, 
cooled  to  room  temperature  and  its  device  characteristics  are  measured.  This 
experiment  will  show  how  the  WSi 2  affects  the  gettering  process. 
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The  other  test  we  are  conducting  checks  the  premise  that  a  number  of  short 
anneals  have  the  same  effect  as  one  long  aggregate.  Towards  that  end  an 
additional  wafer  is  taken,  it  is  given  a  backside  getter,  and  it  is  then 
contacted  with  WSi 2 *  It  is  then  annealed  8  separate  times  for  1  hour  at  a 
time.  At  the  end  of  this  series  of  anneals,  the  wafer's  device  characteristics 
are  compared  to  the  other  two  devices  described  above.  This  experiment  essen¬ 
tially  determines  whether  or  not  the  getter ing  process  can  be  treated  as  a 
linear  system.  If  it  can,  then  the  device  characteristics  after  8  one  hour 
anneals  will  be  indistinguishable  from  one  8  hour  anneal. 

In  summary,  we  have  shown  that  our  measurements  of  device  characteristics 
do  correlate  with  levels  of  metallic  contamination.  The  device  characteristics 
seem  to  be  unaffected  by  vacancies  in  silicon,  and  because  of  that  fact  the 
device  characteristics  are  relatively  unaffected  by  such  parameters  as  pull  rate 
and  cooling  rate  from  high  temperatue  anneals.  With  that  assurance,  we  went 
ahead  and  have  fabricated  a  large  number  of  wafers  with  various  test  structures. 

Using  these  wafers  a  number  of  experiments  have  been  started  checking  some 
of  the  underlying  assumptions  in  our  experimental  procedure. 
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temperature  T. 
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10.  Si  WAFER  SURFACE  PROPERTIES:  SURFACE  COMPOSITION  AND  THE  EFFECT  OF 
CLEANING  PROCEDURE  STUDIED  BY  ELLIPSOMETRY  AND  AUGER  ELECTPON  SPECTROSCOPY 

M.  TAUBENBLATT  and  C.  HELMS 

In  this  program,  we  have  begun  an  investigation  of  the  chemistry  of  wafer 
surfaces  and  the  effect  of  various  cleaning  steps  on  the  chemistry.  One  goal  of 
this  study  is  to  compare  ellipsometry  measurements  to  electron  spectroscopy  to 
provide  more  quantitative  measurements  for  thin  oxides. 

Although  ellipsometry  does  not  provide  an  actual  identification  of  the 
species  present  on  the  surface,  it  does  give  information  concerning  the 
thickness  and  index  of  refraction  of  surface  layers  present  (if  the  optical 
properties  of  the  substrate  are  known).  One  difficulty  with  the  technique  is 
that  the  "effective  substrate  index"  is  a  strong  function  of  surface  damage  and 
adsorbed  layers  such  as  H20,  OH,  etc.,  so  that  care  must  be  taken  to  assure  the 
value  of  substrate  index  used  is  appropriate.  The  procedure  normally  followed 
in  most  laboratories  is  to  fix  the  Si  index  at  3.85  -  0.02  i  and  using  measured 
values  of  a  and  4,  vary  the  layer  thickness  and  index  until  the’ best  fit  to  the 
ellipsometry  equations  is  obtained.  If  the  substrate  index  is  correct,  this 
leads  to  relatively  good  values  for  the  layer  thickness  for  thin  layers  but 
rather  uncertain  values  for  the  layer  index  due  to  the  insensitivity  of  a  and  4, 
to  variations  in  the  index.  It  is  also  important  to  note  that  is  relatively 
insensitive  to  thickness  for  thin  layers,  whereas  a  is  nearly  linear;  a  is 
therefore  much  more  reliable  for  thickness  determination.  Also,  a  and  exchange 
very  little  with  the  layer  index  for  thin  oxides  so  that  determinations  of  the 
layer  index  from  measurements  of  a  and  4,  is  very  inaccurate.  The  importance  of 
the  substrate  index  can  be  seen  if  we  look  at  a  change  in  the  imaginary  part  of 
the  substrate  index  from  -0.02  to  -0.2;  this  gives  a  change  in  a  of  -  7°. 
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A  7°  error  in  a,  however,  gives  a  25  A  error  in  thickness  (a  large  error, 
indeed,  for  thin  oxides).  Our  approach  is,  therefore,  to  rely  on  values  of  a 
starting  from  a  measurement  of  the  substrate  after  removal  from  HF  and  then 
compare  to  the  value  obtained  from  the  wafer  with  the  native  oxide  present. 
Average  values  for  o  averaged  over  a  number  of  samples  for  the  substrate  after 
an  HF  dip  are  «  176.4°  and  after  an  HF  dip  plus  3  hours  are  175.9  [<•  2  A).  As 
received  were  measured  within  30  minutes  of  removal  from  the  shipping  container 
and  give  values  of  a  173.8  (9  A).  The  substrat  value  for  a  of  176.4°  corres¬ 
ponds  to  an  imaginary  part  of  the  substrate  index  of  -0.08  or,  conversely,  7  A 
of  native  oxide  which,  from  our  electron  spectroscopy  studies,  can  certainly  be 
shown  not  to  be  present. 

Unlike  ell ipsometry.  Auger  electron  spectroscopy  provides  detailed  infor¬ 
mation  on  the  constituents  present  within  5-50  A  of  the  surface.  In  addition, 
information  concerning  the  bonding  of  the  surface  atoms  is  also  accessible. 
Measurements  on  wafers  as  received  are,  except  for  the  oxides  present,  quite 
clean  with  the  average  carbon  level  observed  being  0.10  monolayers  in  the 
surface  region.  The  oxygen  concentration  is  approximately  2  monolayers  (5  A). 

In  comparing  the  ellipsometry  results  to  the  Auger  electron  spectroscopy 
results,  we  see  first  that  the  layer  thickness  obtained  from  ellipsometry  is 
larger  («  a  factor  of  two)  than  that  obtained  by  AES.  This  discrepancy  could  be 
caused  by  a  number  of  factors.  First,  due  to  the  uncertainty  in  the  effective 
substrate  index,  a  large  absolute  error  in  the  ellipsometry  measurements  might 
be  expected.  If  the  surface  oxide  is  homogenous,  the  10  A  +  values  of  surface 
oxide  thickness  are  definitely  inconsistent  with  the  AES  results.  It  is 
possible,  however,  that  the  surface  oxide  is  nonuniform  so  that  there  are 
substantial  regions  of  thinner  oxide  present  on  the  surface.  This  could  account 
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for  some  of  the  discrepancy  as  ellipsometry  and  AES  are  effected  much  differ¬ 
ently  by  inhomogenities. 

In  order  to  identify  important  effects  in  the  cleaning  procedures  used 
commonly,  we  subjected  the  wafers  to  various  cleaning  steps  and  then  followed 
the  surface  concentration  with  AES.  A  summary  of  the  results  are  shown  in  Table 
I.  Starting  from  the  bottom,  we  note  that  the  effect  of  a  final  HF  dip  yields 
what  we  believe  to  be  the  undesirable  result  that  the  surface  ends  up  with  a 
very  high  surface  carbon  concentration.  If,  however,  the  HF  dip  is  followed 
(line  5)  by  an  H202  treatment  (note  that  all  these  treatments  were  for  para¬ 
meters  of  the  so-called  "RCA  clean")  without  exposure  to  air  quite  low  carbon 
concentrations  were  obtained.  We  feel  that  the  HF  dip  removes  a  protective 
oxide  layer,  and  when  the  wafer  is  then  exposed  to  air  carbon  from  CO,  C02, 
hydrocarbons,  etc.,  present  can  more  easily  deposit  on  the  surface.  The 
cleanest  surface  was  obtained  using  the  procedure  shown  through  line  5;  however, 
it  may  be  more  desirable  to  leave  the  surface  with  a  thicker  oxide  but  with  a 
minimum  carbon  concentration  (through  line  4  which  is  the  RCA  clean  minus  the 
final  HF  dip). 
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TABLE  I.  EFFECT  OF  VARIOUS  CLEANING  PROCEDURES 


Procedure 


Oxygen 

Concentration 


(1)  As  received 


(2)  H2  SO  4 : H202 


->  2.8  ml 


DI  water 


v 


(3)  NH4OH : H202 : H20 


air 


->  2.3  ml 


V 

DI  water 


(4)  HC1:H202:H20 


V 

DI  water 


air 


->  2.1  ml 


V  air 

(5)  HFDjp  +  DI  -►  H202  -»■  DI  - >  1.1  ml 

|  water  water 

V 

Air 


DI  water 


air 


->  0.4  ml 


Carbon 

Concentration 

0.2  ml 


0.5  ml 

0.1  ml 

0.2  ml 

1.0  ml 


-88- 


11.  STUDIES  OF  OXYGEN  DIFFUSION  THROUGH  Si02 
C.  HELMS 


In  a  recent  paper  Irene  [1]  has  reported  on  measurements  of  the  oxidant 
diffusion  coefficient  using  a  time  lag  method  with  in  situ  ell ipsometry.  For 
the  temperature  range  investigated  (600°  C  -  1000°  C),  a  measurable  lag  time  for 
an  initial  1  pm  thick  oxide  was  observed  only  for  the  1000°  C  sample.  A  value 
of  D  =  2  x  10“13  cm2/sec  was  obtained  which  leads  to  a  value  of  approximately 
1021  cm-3  for  the  oxygen  solubility.  The  value  for  D  is  somewhat  higher  than 
that  obtained  using  exchange  methods  to  measure  diffusion  coefficients  [2-5]  and 
considerably  lower  than  that  obtained  by  Norton  [6]  in  his  permeation  experiment 
and  definitely  larger  than  the  1020  cm-3  maximum  we  have  previously  reported 
[7]. 

As  pointed  out  by  Revesz  and  Schaeffer  [8],  however,  the  comparison  to  the 
exchange  experiments  is  meaningless  in  any  case  since  the  rate  limiting  step  in 
that  process  is  the  exchange  reaction  not  the  diffusion  process.  We  are  then 
left  with  the  large  discrepancy  of  over  104  between  Norton  and  Irene.  The  1021 
cm-3  solubility  is  also  much  larger  than  would  be  expected  for  inert  gas 
molecules  of  roughtly  the  same  size  [8]. 

Recent  results  by  Hamasaki  [9]  seem  to  shed  some  light  on  these  discre¬ 
pancies.  Although  much  of  Hamasaki 's  interpretation  may  be  criticized,  his 
experimental  findings  are  difficult  to  refute.  Namely,  he  finds  that,  when 
sequential  oxidations  are  investigated,  esepcially  when  a  temperature  change  is 
involved,  the  initial  growth  rates  for  the  sequential  oxidation  step  depends 
strongly  on  the  time-temperature  history  of  the  sample  such  that,  for  example, 
"when  the  temperature  of  the  first  oxidation  is  higher  than  that  of  the  second 
oxidation,  a  faster  growth  rate  in  the  initial  state  of  the  second  oxidation 
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is  found.  An  initial  slower  growth  rate  is  also  found  when  the  temperature  of 
the  first  oxidation  is  lower  than  that  of  the  second  oxidation"  [9]. 

The  findings  of  Irene  and  Hamasaki  certainly  indicate  that  the  oxidation  of 
silicon  is  a  much  more  complicated  process  than  we  would  like  to  think.  It 
appears  that  thermal  Si 02  is  not  a  static  quantity,  its  properties  depending  on 
its  temperature  history  and  the  pressure  history  of  soluble  molecules  such  as 
O2,  H20,  CO,  C02,  etc.,  that  may  be  present  in  the  ambient.  If  this  is  the 
case,  the  diffusion  coefficient  would  also  be  expected  to  be  a  function  of  time, 
pressure,  and  temperature,  irrespective  of  any  temperature  dependence  due  to  a 
single  diffusion  activation  energy. 

In  Irene's  case,  the  initial  oxides  were  grown  to  an  initial  thickness  of 
1  pm  at  1000°  C.  An  "inert"  anneal  was  then  performed  for  about  20  hours.  The 
temperature  is  then  changed,  equilibrated,  and  the  ambient  was  switched  to 
oxygen.  The  assumption  was  made  that  the  transport  properties  are  unaffected  by 
either  the  anneal  or  the  temperature  change.  Hamasaki' s  results  certainly  call 
the  second  assumption  into  question. 

Even  the  "inert"  anneal  may  not  be  so  inert,  since  the  structure  of  Si 0 2 
(Si-O-Si  bond  angle  dis*  ;  ution,  microvoid  density,  etc.)  may  well  be  affected 
even  if  the  optical  constants  change  little.  In  addition,  there  is  a  finite 
solubility  of  both  Ar  and  N2  in  SiO 2  which  could  be  responsible  for  blocking 
available  sites  for  02  diffusion. 

It  would  appear  for  most  normal  conditions  where  nonsequential  oxidations 
are  involved  that,  not  only  is  the  oxidant  transport  process  in  steady  state, 
but  the  $i02  itself  reaches  a  steady  state  condition.  This  would  explain  the 
agreement  between  Norton's  values  for  diffusion  coefficient  [5],  the  expected 
solubility  of  02  [8],  and  the  parabolic  rate  constant  for  oxidation  [10]. 
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12.  PHOSPHORUS  DIFFUSION  IN  SILICON 
P.  FAHEY,  R.  DUTTON 

Of  the  four  coirmon  dopants  used  in  Si  device  processing,  P,  As,  B  and  Sb, 
phosphorus  exhibits  by  far  the  most  complex  diffusion  behavior  under  typical 
process  conditions.  Consequently,  it  has  proven  the  most  difficult  to  model. 
Since  the  inception  of  SUPREM,  phosphorus  diffusion  has  been  modeled  using  the 
theory  of  Fair  and  Tsai  [1].  Unfortunately,  model  predictions  do  not  correlate 
well  with  experimental  results  under  all  conditions,  and  the  discrepancies  prove 
worse  under  those  conditions  most  important  to  modern  process  technology,  i.e. , 
low  temperatures  and  shallow  diffusions.  In  addition,  since  the  time  it  was 
first  introduced,  a  number  of  objections  have  been  raised  against  the  physical 
basis  of  the  theory  which  indicate  that  the  present  model  cannot  be  extended, 
even  empirically,  to  handle  important  cases  such  as  phosphorus  diffusion  in 
two-dimensions,  phosphorus  diffusion  from  polysilicon,  and  interaction  of 
phosphorus  with  other  inpurities  present  during  processing.  Despite  these 
shortcomings  the  Fair-Tsai  model  remains  the  only  quantitative  model  of 
phosphorus  diffusion  yet  developed. 

In  this  last  time  period  we  have  thoroughly  reviewed  the  results  of  all 
pertinent  experimental  work  concerning  phosphorus  diffusion  and  have  scrutinized 
all  previously  proposed  theories  in  order  to  develop  a  better  model.  The 
results  of  this  effort  have  been  encouraging  and  have  produced  a  consistent 
picture  of  phosphorus  diffusion  in  Silicon.  At  this  juncture  we  are  able  to 
enumerate  In  detail  the  inadequacies  of  the  present  model  and  to  outline  an 
alternative  model  which  Is  both  theoretically  sound  and  more  consistent  with 
experiment.  The  essence  of  our  model  is  that:  1)  phosphorus  diffusion  occurs 
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through  a  two-stream  mechanism  whereby  phosphorus  atoms  can  migrate  as  either  a 
phosphorus-vacancy  (PV)  pair  or  as  a  phosphorus-interstitialcy  (PI)  pair;  2) 
the  large  fraction  of  electrically  inactive  phosphorus  (•  2/3  total  concen¬ 
tration)  present  under  typical  predeposition  conditions  is  due  to  precipitation, 
and  3)  phosphorus  diffusion  generates  excess  Si  interstitials  which  can  affect 
the  diffusion  of  other  iirpurities.  This  model  differs  markedly  from  the  Fair- 
Tsai  theory.  We  would  like  to  first  justify  the  above  statements  and  contrast 
them  with  the  equivalent  stipulations  of  the  Fair-Tsai  model  in  order  to  eluci¬ 
date  the  differences  in  predictive  capabilities  between  our  model  and  that  of 
Fair  and  Tsai.  We  will  then  discuss  the  experimental  work  which  is  both  in 
progress  and  under  consideration  to  prove  the  correctness  of  our  model. 

We  are  led  to  point  1)  above  by  considering  that  the  kinked,  two-zoned, 
profile  of  phosphorus  diffusion  is  characteristic  of  a  process  where  the 
species  under  consideration  can  exist  in  two  states  and  diffuse  by  two  different 
mechanisms.  The  kinked  profile  is  a  result  of  superposition  of  the  two  diffu¬ 
sive  fluxes.  For  example,  tellerium  in  germanium  (which  exists  in  the  crystal 
in  both  substitutional  and  interstitial  states)  exhibits  a  diffusion  profile 
very  similar  to  that  of  phosphorus.  The  Fair-Tsai  model  proposes  that  phosphorus 
diffuses  primarily  as  P+  V=  pairs.  As  the  local  Fermi  level  (which  depends  on 
the  local  P  concentration)  drops,  the  P+  Vs  pairs  dissociate,  thereby  genera¬ 
ting  excess  vacancies  which  are  responsible  for  formation  of  the  kink.  While 
this  mechanism  is  not  impossible,  it  is  open  to  many  questions  and  will  be 
discussed  shortly.  Our  conclusion  2)  is  based  on  the  direct  observation  by 
Armigliato  et.al.  of  precipitates  [2]  in  concentrations  high  enough  to  account 
for  the  electrically  Inactive  fractions.  Fair  and  Tsai  have  proposed  that  this 
Inactive  fraction  Is  due  entirely  to  (electrically  Inactive)  P+  V*  pairs.  This 
hypothesis  is  shown  here  to  be  Incorrect  on  fundamental  grounds.  Conditions  of 
local  charge  neutrality  dictate  that: 
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n  +  [P+V“]  +  [V]  +  2[V“]  -  p  +  [P+]  +  [P+V°] 


(1) 


while  conservation  of  matter  yields 

(p)totll  ■  [p*l  *  [fV]  *  [pVj.  [pV]  (2) 

(here  brackets  denote  concentrations  and  n  and  p  are  the  electron  and  hole 
concentrations  respectively).  Combining  these  two  equations  we  show  that  for 
the  inactive  concentration 

[P]total  '  n  *  2tp+v=3  +  [p+H  +  2tV=]  +  tV']  >  2[p+V=]  (3> 

Thus,  P+V=  pairs  can  account  for  at  most  one  half  of  the  inactive 
phosphorus.  To  estimate  the  absolute  concentration  of  P+V=  pairs  we  write 
the  formation  reaction  as 


P+  +  V  =  ^  P+  V= 


(4) 


Applying  the  law  of  mass  action  in  a  probabilistic  sense  the  above  reaction 
implies  that  in  equilibrium 


[p+]  [v3] 

[Si]  [Si] 


exp(Efa/kT)  = 


[pV] 

IsiT 


(5) 


where  0  is  a  geometrical  factor  of  order  unity,  [Si]  is  the  concentration  of 
silicon  lattice  sites  (*  5.5  x  1022  cm3),  and  E&  is  the  binding  energy  of  the 
P+Va  pair.  Using  the  commonly  accepted  vacancy  energy  level  of  E3  *  Ec 
-  O.lleV  for  the  double  negatively  charged  vacancy,  and  the  highest 
experimentally  observed  value  of  n  (»  4  x  1020  cnr3),  we  obtain  for  1000°  C 

[pV] 

-  *  3  X  1010  exp(Eh/kT)  (6) 

[P+] 
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Under  all 


For  Eb  as  high  as  l.OeV  we  obtain  [P+V*]/[p+]  -  5  x  10"5  cm-3, 
realizable  conditions  then  [P+V“]  «  [P+].  In  terms  of  absolute 
concentration  [P+V®]  is  probably  in  the  range  of  10 13  to  10 16  cm"3 

Our  claim  3)  that  phosphorus  diffusion  generates  excess  Interstitials  is 
supported  by  the  facts  that,  the  diffusivityof  a  boron  base  is  considerably 
enhanced  (emitter  dip  effect),  and  surface  stacking  faults  (which  are  extra 
planes  of  atoms),  grow  during  the  diffusion  of  phosphorus.  It  is  generally 
acknowledged  that  these  same  effects  occur  during  oxidation  of  the  silicon 
surface  due  to  the  generation  of  excess  interstititals.  Our  proposition  is  in 
direct  opposition  to  the  Fair-Tsai  theory  which  proposes  the  generation  of 
vacancies  during  diffusion  (in  this  case  one  would  expect  enhanced  shrinkage 
rather  than  growth  of  stacking  faults).  We  also  note  that  since  the  phosphorus 
precipitates  observed  by  TEM  experiments  are  of  an  interstitial  nature,  it  is 
natural  to  expect  the  precipitation  process  itself  to  generate  excess  Si 
interstitials.  This  is  in  accord  with  experimental  results  of  diffusion  of 
phosphorus  from  polysilicon  in  which  precipitation  occurring  primarily  in  the 
polysilicon  is  coincident  with  a  reduction  of  the  emitter  dip  effect.  This  type 
of  effect  is  impossible  to  model  with  the  Fair-Tsai  theory. 

To  verify  our  model  exper imentally  we  have  attempted  to  observe  the 
behavior  of  antimony  buried  layers  during  the  in-diffusion  of  phosphorus  at  the 
surface.  This  experiment  is  based  on  the  following  considerations.  Due  to  its 
size,  i.e.,  greater  than  Si,  antimony  is  expected  to  have  a  relatively  small 
interstitialcy  component  of  diffusion  and  to  effect  movement  primarily  through  a 
vacancy  mechanism.  Consistent  with  this  idea  is  the  experimental  observation 
that  when  the  silicon  surface  is  oxidized  antimony  diffusion  is  not  enhanced  and 
in  some  cases  It  Is  actually  retarded  (it  is  postulated)  due  to  annihilation  of 
vacancies  by  the  excess  interstitials  generated.  It  is  clear  then  that  antimony 
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should  be  enhanced  during  phosphorus  diffusion  if  vacancies  are  generated,  ana 
should  display  no  enhancement  or  possibly  retardation  if  interstitials  are 
generated.  Preliminary  results  indicate  that  antimony  diffusion  is  indeed 
retarded. 

A  number  of  additional  experiments  are  presently  being  developed  in  an 
effort  to  further  test  the  correctness,  and  quantify,  our  new  model.  It  is 
hoped  that  antimony  can  be  used  under  a  variety  of  conditions  as  a  marker  layer 
to  differentiate  between  vacancy  and  interstitial cy  mechanisms  of  diffusion. 
Experiments  to  detect  the  phosphorus  interstitialcy  directly  (e.g.,  by 
channeling)  are  in  a  preliminary  state,  and  have  thus  far  been  restricted  to 
feasibility  studies. 
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